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ABSTRACT 


The  radar  cross  section  of  a  continuous,  convex,  body  of  revo¬ 
lution  composed  of  N  sections,  each  section  described  by  a  second- 
degree  equation  has  been  analyzed  using  the  geometrical  theory  of  dif¬ 
fraction.  Wedge  diffraction  has  been  applied  to  determine  the  scattered 
field  due  to  discontinuities  in  slope  between  sections  of  the  target,  and 
creeping  wave  theory  has  been  applied  to  determine  the  scattered  field 
due  to  propagation  of  energy  around  the  target.  A  solution  for  the 
diffracted  field  on  an  axial  caustic  is  presented.  An  approximate 
solution  for  the  scattered  field  near  and  at  the  normal  direction  to  a 
conical  generator  is  developed.  A  Simplified  ray  path  geometry11 
for  the  creeping  wave  is  presented  and  related  to  the  scattering  by 
spheres  and  prolate  spheroids.  The  H-plane  field  of  the  sphere  is 
calculated  using  creeping  wave  techniques  for  a  ray  geometry  defined 
by  the  Poynting  vector  at  the  shadow  boundary.  The  approximate 
creeping  wave  solution  for  the  edge-on  backscattering  of  disks  is 
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CHAPTER  I 
INTRODUCTION 


A.  Components  of  the  Backscattered 
Field  from  a  Composite  Target 

The  determination  of  the  backscattered  fields  of  a  composite 
conducting  target  is  a  practical  problem  of  interest.  Such  a  target 
which  may  be  composed  of  smooth  surfaces,  flat  facets,  edges,  fins, 
and  cavities  presents  a  formidable  problem  in  analysis.  An  exact 
solution  using  analytical  techniques  such  as  separation  of  variables 
is,  in  general,  impractical  for  such  a  target.  Thus  the  usual 
approach  to  the  analysis  is  by  summation  of  solutions  for  the  back- 
scatter  of  the  component  parts  of  the  target.  This  approach  has 
enjoyed  good  success  in  determining  the  major  scattering  mechanisms 
for  complex  targets.  However,  previous  analysis  has  been  devoted 
to  specified  targets.  An  approach  which  will  apply  to  a  class  of 
targets  whose  individual  characteristics  may  be  varied  at  will  is 
needed.  Several  techniques  have  been  developed  for  the  analysis 
of  the  general  target.  These  consist  of  the  point-matching  computer 
solution,  and  the  wire  grid  model  computer  solution.  Both  of  these 
approaches  yield  good  results  for  targets  whose  electrical  size  is 
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small,  resulting  in  a  set  of  N  linear  equations  which  are  within  the 
capacity  of  present  digital  computers  to  handle.  The  "Third  Gener¬ 
ation"  computers  with  their  greatly  expanded  capabilities  promise 
extension  of  these  techniques  to  larger  targets.  However,  the  time 
involved  in  such  computations  is  great  resulting  in  a  high  cost  per 
data  point.  Thus  a  relatively  rapid  and  inexpensive  technique  would 
be  valuable,  especially  where  real  time  simulation  of  radar  problems 
is  desirable.  An  analysis  based  upon  geometrical  optics,  and  the 
geometrical  theory  of  diffraction,  is  presented  herein  which  yields 
approximate  results  with  relatively  little  expenditure  of  computer 
time.  This  technique  is  suited  to  real  time  computation  for  simu¬ 
lation  of  radar  problems. 

Thus  the  task  of  this  study  consists  of  assembling  known  solu¬ 
tions  for  the  particular  scattering  mechanisms  involved,  extending 
these  solutions  where  necessary  and  applying  the  solutions  in  the 
development  of  a  computer  program  to  determine  the  monostatic 
radar  cross  section  for  a  large  class  of  axially  symmetric  targets. 
Future  problems  involved  in  extending  the  computer  program  include 
nonsymmetric  targets  and  bistatic  scattering.  These  problems  are 
discussed  briefly  and  it  is  seen  that  such  extensions  are  possible 
using  the  solutions  which  have  been  applied  to  the  restricted  class 
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of  targets.  The  objective  of  this  study  is  to  demonstrate  the  capabili¬ 
ty  of  these  techniques  in  a  practical  manner. 

B.  Description  of  the  Principal 
Backscattering  Mechanisms 

The  principal  backscattering  mechanisms  of  a  composite  target 
can  best  be  illustrated  pictorially.  Figure  1  illustrates  a  composite 
target  viewed  by  an  interrogating  radar  signal.  The  numbered  po¬ 
sitions  on  the  target  correspond  to  different  scattering  mechanisms 


Fig.  I  -  -A  general  target. 
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on  the  target  which  may  contribute  to  the  backscattered  field.  We 
divide  these  mechanisms  into  five  classes  as  follows: 

1.  The  geometrical  optics  field 

The  geometrical  optics  field  results  from  reflection  of  the 
incident  energy  at  the  specular  point  on  the  target.  The  resulting 
contribution  to  the  backscattered  field  may  be  evaluated  using  geo¬ 
metrical  optics  techniques  involving  the  Gaussian  curvature  at  the 
specular  point.  The  specular  point  is  defined  as  that  point  in  the 
illuminated  region  where  the  surface  normal  is  parallel  with  the 
direction  of  incident  wave  propagation. 

The  geometrical  optics  field  solution  fails  for  some  targets 
such  as  the  cone,  cone  sphere  junction,  and  flat  plates.  For  such 
targets  the  physical  optics  solution  is  required. 

2.  Wedge  diffracted  fields 

The  wedge  diffracted  fields  result  from  diffraction  at  slope 
discontinuities  on  the  target.  This  contribution  can  be  evaluated 
using  the  techniques  of  the  Geometrical  Theory  of  Diffraction 
developed  by  Keller[  l]  . 

The  effects  of  fins  at  aspects  removed  from  the  normal  to  the 
plane  of  the  fin  may  be  evaluated  using  wedge  diffraction  techniques. 


Normal  to  the  plane  of  the  fin  the  techniques  of  physical  optics  may 
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applied.  However,  the  vertices  of  fins  cannot  be  treated  using  the  geo¬ 
metrical  theory  of  diffraction  as  the  diffraction  coefficient  of  a  vertex 
is  not  known. 

3.  Creeping  wave  fields 

The  creeping  wave  backscattered  fields  arise  from  that  portion 
of  the  incident  energy  which  is  trapped  at  the  surface  of  the  target, 
and  propagates  around  the  target,  eventually  reradiating  energy  in  the 
backscatter  direction.  These  fields  may  be  computed  using  the  geo¬ 
metrical  theory  of  diffraction  and  knowledge  of  the  differential  geome¬ 
try  of  the  target  surface. 

4.  Cavity  fields 

The  cavity  fields  are  caused  by  concavities  in  the  target  which 
may  focus  energy  as  for  a  corner  reflector.  The  point  matching  tech- 

nique[  Z]  has  been  applied  to  the  determination  of  cavity  contributed 
fields.  In  some  cases  these  cavity  contributions  can  be  analyzed  as 
antenna  contributions. 

5.  Antenna  mode  fields 

The  antenna  contributions  are  due  to  antennas  on  the  target  which 
receive  and  then  reradiate  energy  in  the  backscatter  direction.  These 
contributions  can  be  treated  using  antenna  scattering  technique  s[  3 ,  4]  . 


6 


It  is  apparent  that  an  analysis  which  combines  all  the  scattering 
mechanisms  listed  above  would  be  a  truly  formidable  task.  It  is  the 
purpose  of  this  study  to  treat  the  contributions  due  to  the  first  four 
mechanisms,  i.  e.  ,  geometrical  optics  fields,  wedge  diffracted  fields, 
fin  diffracted  fields,  and  creeping  wave  fields,  for  a  general  convex 
second  degree  surface  of  revolution. 

C.  Superposition  of  the  Geometrical 
Optics,  Wedge  Diffraction  and 

Creeping-Wave  Contributions 

In  order  to  obtain  a  solution  for  the  backscattered  fields,  the 
contributing  mechanisms  and  their  locations  on  the  target  must  be 
identified.  Thus  it  is  necessary  to  scan  the  surface  of  the  target  to 
obtain  the  specular  point,  the  location  of  wedges,  and  the  points  of 
attachment  and  reradiation  of  the  creeping  waves.  This  portion  of 
the  solution  thus  deals  with  the  geometrical  properties  of  the  target. 

In  particular,  the  normal  vector  to  the  surface  must  be  calculated  in 
order  to  determine  the  location  of  the  specular  point  and  the  attach¬ 
ment  and  reradiation  points  of  the  creeping  waves.  The  normal 
vector  is  also  needed  to  identify  the  location  and  included  angle  of 
the  wedges.  In  order  to  calculate  the  creeping  wave  path  lengths 
the  differential  arc  length  on  the  surface  is  needed  as  well  as  the 
principal  radii  of  curvature  at  each  point  along  the  path.  The 
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Gaussian  curvature  at  the  specular  point  is  needed  to  calculate  the 
geometrical  optics  field. 

Given  the  geometrical  properties  of  the  surface  of  the  target, 
the  solutions  for  the  scattered  fields  must  next  be  obtained.  The 
specular  scattered  field  is  easily  obtained  from  the  Gaussian 
curvature  at  the  specular  point.  The  wedge  diffracted  fields  may 
be  obtained  using  the  Keller[  l]  or  Pauli[  5]  solution  for  we  dge 
diffraction  and  the  ray  techniques  of  the  Geometrical  Theory  of  Dif¬ 
fraction.  However  the  creeping  wave  contributions  are  not  easy  to 
obtain  since  the  complete  creeping  wave  solution  is  available  only 
for  the  cylinder  and  sphere.  The  behavior  of  the  creeping  wave  dif¬ 
fraction  and  attenuation  coefficients  and  the  geodesic  ray  paths  are 
not  known  for  a  general  body.  Thus  it  is  necessary  to  obtain  an 
approximate  solution  for  the  behavior  of  a  creeping  wave  on  a  general 
surface.  The  effects  of  ray  path  geometry  on  the  sphere  have  been 
examined  in  order  to  obtain  an  approximate  simplified  ray  path  for 
a  general  target.  Diffraction  and  attenuation  coefficients  have  been 
developed  empirically  for  such  a  simplified  ray  path  and  have  been 
shown  to  yield  good  results  for  prolate  spheroids,  ogives,  and  disks. 
The  assumption  of  a  simplified  ray  path  on  the  surface  reduces  the 
computation  complexity  and  allows  one  to  reduce  the  determination 
of  the  ray  path  on  a  general  target  to  a  single  numerical  iteration  of 
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the  differential  equation  for  the  geodesic  on  a  surface.  In  addition, 
such  a  simplified  ray  picture  allows  one  to  take  into  account  the 
effects  of  intersections  on  the  surface  which  have  discontinuities  in 
the  first  derivative  along  the  ray  path.  These  effects  are  treated 
only  in  the  first  order  in  this  study  but  extension  to  higher  order 
effects  is  straightforward. 

The  determination  of  the  backscatter ed  fields  from  a  general 
target  is  thus  accomplished  by  utilizing  geometrical  optics  and  wedge 
diffraction  techniques  for  the  illuminated  portions  of  the  target  and  an 
approximate  creeping  wave  solution  for  the  contributions  from  the 
shadow  zone  of  the  target.  A  computer  program  coded  in  Fortran  IV 
is  the  end  result  of  this  investigation.  This  program  computes  the  back- 
scattered  fields  of  a  general  second  order  surface  of  revolution  using 
geometrical  optics,  wedge  diffraction,  and  creeping  wave  techniques 
for  the  case  of  parallel  polarization,  and  geometrical  optics  and 
wedge  diffraction  techniques  for  the  case  of  perpendicular  polari¬ 
zation  of  the  incident  field.  The  program  is  suited  for  on-line  simu¬ 
lation  of  radar  problems.  This  program  has  been  tested  for  both 
canonical  shapes  and  complex  targets  and  the  results  have  been  com¬ 
pared  to  measured  data  of  test  targets.  Results  to  date  indicate  that 
the  computer  solution  is  accurate  to  within  3  dB  of  the  measured 
data  in  the  regions  where  the  solutions  applied  are  valid. 


CHAPTER  II 
GEOMETRICAL  OPTICS 


A.  Ray  Optics 

The  analysis  of  wave  propagation  using  geometrical  optics  is 
based  upon  the  assumption  that  energy  travels  along  straight  lines, 
called  rays,  except  when  modified  by  reflection  or  refraction.  Or¬ 
thogonal  to  the  ray  trajectories  there  exists  a  set  of  equiphase 
surfaces.  Figure  2  illustrates  this  orthogonal  coordinate  set.  The 
relation  of  geometrical  optics  to  electromagnetic  theory  has  been 
derived  by  Luneberg[  6]  through  the  application  of  Maxwell's  e- 
quations  and  the  boundary  conditions.  For  the  purposes  of  this 
study  a  brief  description  of  the  methods  required  to  calculate  the 
scattered  fields  of  interest  will  be  presented  here. 

In  an  isotropic,  homogeneous  medium  the  rays  are  straight 
lines,  and  the  law  of  reflection  and  Snell's  law  of  refraction  de¬ 
scribe  the  behavior  of  the  rays  at  a  boundary. 

(  1)  0O  =  01  (  Law  of  Reflection) 

(2)  'Ju  sin  0Q  =  sr2  sin  G2  (Snell's  Law), 

where  ,  e2  are  the  permittivities  of  the  media. 
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Fig.  2--The  relation  of  rays  and  wavefronts. 

Figure  3  illustrates  these  laws;  the  incident,  reflected,  and  re¬ 
fracted  rays  and  the  normal  to  the  surface  at  the  point  of  encounter 
are  coplanar.  In  the  investigation  of  scattering  by  conducting  bodies 
the  law  of  reflection  is  used.  However  scattering  by  penetrable 
bodies  such  as  dielectric  spheres  requires  use  of  Snell’s  Law. 

Such  analyses  have  been  performed  by  Peters  and  Thomas[  7]  , 
Kouyoumjian,  Peters,  and  Thomas[  8]  ,  Swarner  and  Peters[  9]  , 
and  Peters,  Kawano,  and  Swarner[  1 0 ]  . 
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Fig.  3- - Reflection  and  refraction  at  a  boundary 
between  two  media 

In  addition  to  specifying  the  ray  trajectories  it  is  necessary  to 
account  for  the  amplitude  behavior  of  the  field.  This  is  accomplished 
by  applying  conservation  of  energy  within  the  astigmatic  flux  tube 
depicted  in  Fig.  4.  It  is  assumed  that  the  field  associated  with  a 
particular  ray  u(  l  )  can  be  described  as 

(3)  u(i  )  =  A0ej4>  F(i  )  e"jki 

whe  re 

A 0  is  the  amplitude  at  a  reference  point, 

4>  is  the  phase  at  the  reference  point, 

F(  f  )  is  the  spatial  attenuation  factor,  and 

e-jW 


is  the  phase  factor  . 
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Fig.  4- -An  astigmatic  ray  tube 


F(  !L  )  is  dependent  upon  the  nature  of  the  reference  equiphase  surface 
and  accounts  for  the  convergence  or  divergence  of  the  flux  tube  in  the 
direction  H  .  The  flux  tube  depicted  in  Fig.  4  is  used  to  determine 
F(  i  )  as  follows.  It  is  seen  that  the  principal  radii  of  curvature  of  the 
cross  section  dcrQ  are  p*  and  p2.  The  field  amplitude  at  the  reference 
point  O  is  taken  to  be  A0?  and  the  field  amplitude  at  the  distance  H  is 
A.  The  energy  in  the  wave  is  proportional  to  the  square  of  the 
amplitude.  Applying  conservation  of  energy  within  the  flux  tube  we 
have 

2  2 

(  4)  A0dcr0  =  A  dcr  =  constant 

where  dcr  is  the  c  ross  section  of  the  flux  tube  at  i  ,  having  principal 
radii  of  curvature  px  +  l  and  p2  +  H  . 
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The  ratio  of  their  areas  can  be  written  as 


(5) 


do-  _  (  px  +  l  )  (  P2  +  i  ) 


dcr 


o 


Pi  P2 


Inserting  this  result  into  Eq.  (  4)  yields 


A 


Pi  P2 


(  6) 


\  (  Pi  +  *  )  (  P2 +  i  ) 


The  field  u(  i  )  is  then 


(V 


u(  i  )  =  A0e^° 


Pi  P2 


e 


-jki 


\  (  Pl  +  ^  )  (  P2+  *  ) 


At  the  locations  1  -  -  pi  and  i  =  -p2  the  field  given  by  Eq.  (  7)  becomes 
infinite,  and  the  ray  optics  solution  fails  to  obtain  the  correct  value 
of  u(  -pi)  or  u(  -p2).  These  locations  are  termed  "caustics"  of  the 
geometrical  field.  In  order  to  calculate  the  field  at  such  a  caustic 
point  Kay  and  Kellerf  1  1  ]  have  derived  a  caustic  correction  factor. 
The  derivation  of  such  a  correction  factor  proceeds  from  a  solution 
of  the  scalar  wave  equation,  and  will  not  be  presented  here.  How¬ 
ever,  certain  results  of  this  derivation  are  necessary.  It  has  been 
demonstrated  that  in  the  region  away  from  the  caustics  the  geo¬ 


metrical  optics  solution  is  correct  if  a  phase  shift  of 


introduced  upon  traversal  of  a  caustic  line.  This  phase  shift  is 
accounted  for  by  Eq.  (  7)  if  the  sign  of  I  is  preserved.  If 
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-  p2  <  i  <  -pj  ,  i.  e.  ,  the  caustic  line  has  been  crossed,  then 


Pl  +  1  <  0,  p2  +  i  >  0,  and  u(  1  )  is 


(  8) 


Pl  P2 


e 


-jki 


\  -(  Pl  +l  )  (  P2  +  i  ) 


-iz  e-jkf 


TT 


Pl  P2 


e 


^  (  Pl  +  t  )  (  P2+  t  ) 


and  the  phase  shift  is  apparent. 

If  a  ray  strikes  a  boundary,  and  is  transformed  into  reflected 


and  refracted  rays  with  directions  specified  by  the  laws  of  reflection 


and  refraction  the  values  of  A0,  Pi  »  and  p2,  must  be  determined 

for  both  the  reflected  and  refracted  rays.  In  the  general  case  the 
reflection  and  transmission  coefficients  are  polarization  sensitive 
except  for  normal  incidence.  This  polarization  sensitivity  does  not 
concern  us  here  as  the  targets  to  be  treated  are  conducting  bodies. 
However,  in  the  case  of  penetrable  bodies  these  coefficients  must  be 
evaluated  (  see  for  example  Thomas[  1 Z]  ). 

B.  Specular  Scattering  of  the 
Reflecting  Surface 

Consider  a  ray  incident  upon  a  curved  surface  as  shown  in 
Fig.  5.  The  two  dimensional  case  is  illustrated  for  simplicity. 

The  source  of  cylindrical  rays  is  a  point  P  located  a  distance  S.  Q 


from  point  Q,  the  intercept  of  the  rays  and  the  reflecting  surface  S. 
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Fig.  5- -R eflection  by  a  singly  curved  surface 


The  surface  has  a  radius  of  curvature  rs  at  point  Q,  and  the  angle  be¬ 
tween  the  incident  ray  and  the  surface  normal  at  Q  is  0Q.  It  is  now  a 
geometrical  problem  to  determine  using  the  law  of  reflection,  the 
direction  and  divergence  of  the  reflected  rays.  In  order  to  determine 
the  spatial  attenuation  factor  two  rays  originating  at  P  are  used. 

These  rays  have  a  small  angular  deviation  (3.  Provided  the  angular 
difference  Ao  of  the  reflected  rays  is  small  we  have 


A(0o-a) 
cos  0Q 


A(0o  +  a) 

-  =  rs  * <* 

o 


(9) 


=  Pi 


cos  0 
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Upon  reduction  one  has 


(  10) 


JL^ 

Pi 


Thus  we  have  obtained  the  caustic  distance  pj  of  the  virtual  focus  in 
terms  of  the  local  radius  of  curvature  of  the  surface,  the  reflection 
angle  and  the  caustic  distance  of  the  incident  ray. 

A  similar  situation  exists  for  the  case  of  a  three-dimensional 
problem.  A  three-dimensional  problem  will  involve  two  caustics  of 
the  reflected  ray  tube  (which  may  not  be  coincident)  similar  to  the 
case  shown  in  Fig.  6.  The  same  procedure  may  be  used  to  obtain 
the  caustic  distances  by  separating  the  problem  into  two  cylindrical 
problems.  In  this  case  the  two  principal  radii  of  curvature  of  the 
surface  at  the  reflection  point  are  required.  The  principal  radii  of 
curvature  may  be  used  in  conjunction  with  Euler's  theorem  to  de¬ 
termine  the  radius  of  curvature  in  a  direction  which  is  not  coincident 
with  a  coordinate  direction  as  noted  in  Appendix  L 

Once  the  reflection  point,  angle  of  reflection  and  the  local 
radius  of  curvature  are  known  the  divergence  factor  can  be  de¬ 
termined.  This  information  in  connection  with  Eq.  (  10)  allows  the 
reflected  field  to  be  calculated.  In  a  later  chapter  it  will  be  demon¬ 
strated  that  these  techniques  are  more  generally  applicable,  specifi¬ 
cally  to  the  case  of  diffraction  by  a  curved  edge  in  connection  with 


the  geometrical  theory  of  diffraction. 
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Fig.  6-~Reflection  by  a  double- cur ved  surface* 

Consider  the  situation  depicted  in  Fig.  6.  A  spherical  bundle  of 
rays  originating  at  P  is  incident  upon  a  three-dimensional  surface  S 
at  the  reflection  point  Q.  It  is  desired  to  obtain  the  scattered  field 
through  the  use  of  ray  optics  techniques.  Let  the  coordinate  system 
describing  the  surface  be  the  spherical  system  (  r,  0,4>).  The  rays 
bounding  the  incident  flux  tube  can  be  taken  to  lie  in  the  0  and  4> 
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planes  respectively  with  no  loss  in  generality.  Thus  the  incremental 
area  of  the  surface  about  the  reflection  point  can  be  written  as 


(  11)  do  Q  =  r2sin  0  d0  d<j> 


We  may  now  separate  the  problem  into  two  two-dimensional  problems. 

Referring  to  Eq.  (  10)  we  may  determine  the  caustic  distances  in 

G  cb 

terms  of  the  principal  radii  of  curvature  r  ,  rY  of  the  surface  as 


(  12) 


1  Z  1 

—  — £  + 

pi  r  cos  0O  iQ 


(  13) 


1 _ +  J_ 


P2  r^cos  4>q  iQ 
where  0O  and  §Q  are  defined  as  in  Fig.  5  (  i.  e.  ,  they  are  not  polar 

angles)  .  Thus  we  may  write  the  reflected  (  i.  e.  ,  scattered)  field  as 

'  ^/^ 


(14)  u(  i)  =  A0e 


j^o 


Pi  P2 


(  Pi  +  *  )  (  P2+  *  ) 


-jki 


The  case  of  backscatter  (  i.  e.  ,  GQ  -  cj>0  -  0)  results  in 


(  15) 


u(  &  )  =A0e 


j^o 


r^rVn 


(  r0io+  l  )  (  r4>i0+  i  )  (  2i0+  r^)  (  2i0  +  r6) 


1/2 


-jki 


For  the  case  of  an  incident  plane  wave  (i  — *  °°)  the  far-zone  scat¬ 


tered  field  (  i  >  >  0)  is 
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(  16) 


e 


l 


-jki 


The  spatial  attenuation  factor  appearing  in  Eq.  (15)  is  seen  to  be  one 
half  the  reciprocal  of  the  square  root  of  the  Gaussian  curvature  of  the 
surface  at  the  reflection  point.  Thus  the  problem  of  determining  the 
geometrical  optics  backscatter  from  a  three-dimensional  target  re¬ 
duces  to  the  differential  geometry  problem  of  the  determination  of 
the  Gaussian  curvature  given  in  Appendix  L 


CHAPTER  IE 
WEDGE  DIFFRACTION 


A.  Single  Diffraction 

Consider  the  case  of  a  perfectly-conducting  wedge  illuminated  by 
a  monochromatic  plane  wave  as  illustrated  in  Fig.  7.  According  to 


Fig.  7--Geometrical  optics  rays  for  cylindircal 
wave  incidence  on  a  wedge 

the  principles  of  geometrical  optics  two  classes  of  rays  may  exist. 
These  are  the  direct  rays  A-B  and  the  reflected  rays  A-C-B.  The 
behavior  of  these  rays  is  determined  by  Fermat's  principle  which 
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may  be  stated  as  follows:  MThe  time  elapsed  in  the  passage  of  light 
between  two  fixed  points  is  an  extremum  with  respect  to  possible 
paths  connecting  the  points.  n  Or  equivalently  that  the  value  of  the 
integral  between  points  Pi  ,  P2  on  a  path  s 


where  n  is  the  index  of  refraction  of  the  medium,  be  an  extremum. 
Keller[  1]  has  extended  Fermat’s  principle  to  edge  diffraction  in 
three  dimensions  through  the  following  assumption:  nA  singly  dif¬ 
fracted  ray  connecting  two  points  is  a  curve  whose  length  is  station¬ 
ary  among  all  curves  connecting  these  two  points  and  having  one 
point  on  the  edge.  M  Asa  consequence  of  this  extended  principle, 
a  ray  normally  incident  upon  the  edge  of  a  wedge  generates  a  family 
of  diffracted  rays  which  lie  in  a  disc  having  the  edge  as  its  axis. 
Also  a  ray  incident  obliquely  to  the  edge  of  a  wedge  generates  a 
family  of  diffracted  rays  which  lie  on  a  cone  having  the  edge  as  an 
axis  and  a  half-angle  of  the  cone  equal  to  the  angle  between  the 
incident  ray  and  the  edge.  These  two  cases  are  shown  in  Figs.  8 
and  9.  In  addition  Kellerf  l]  has  shown  that  for  diffraction  in  three 
dimensions  the  two-dimensional  diffraction  coefficients  must  be 
modified  to  take  into  account  the  distribution  of  energy  in  the  cone 


of  diffracted  rays. 
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Diffraction  of  a  cylindrical  wave  by  a  wedge. 


WEDGE 


INCIDENT 

RAY 


7>- 


EDGE 


DIFFRACTED 

RAYS 


Fig.  9- - Three- dimensional  picture  of  the  rays 
diffracted  by  a  wedge. 


24 


In  the  geometrical  theory  of  diffraction  the  use  of  ray  techniques 
to  describe  energy  flow  is  employed  in  the  same  manner  as  in  geo¬ 
metrical  optics.  In  fact,  the  treatment  is  the  same  as  in  geometrical 
optics,  except  for  diffraction  effects  which  are  caused  by  discontinuities 
of  the  structures  involved.  In  conventional  geometrical  optics,  illumi¬ 
nation  of  structural  discontinuities  causes  shadow  boundaries  about 
which  the  field  is  discontinuous;  its  unperturbed  value  is  on  the  illumi¬ 
nated  size  of  the  shadow  boundary  and  the  null  value  is  on  the  shadow 
side.  A  similar  situation  exists  for  the  reflected  rays  for  which  there 
is  also  a  shadow  boundary.  The  geometrical  theory  of  diffraction  takes 
into  account  the  diffraction  effect  of  a  structural  discontinuity  by  using 
solutions  of  canonical  problems,  and  expresses  the  field  more  exactly 
by  eliminating  the  apparent  discontinuity  of  the  field  at  the  shadow 
boundaries  introduced  by  geometrical  optics.  In  particular,  the  dif¬ 
fraction  by  a  wedge  of  perfect  conductivity  is  one  of  the  most  important 
canonical  problems  and  is  the  one  employed  in  this  research.  The  dif¬ 
fracted  field  introduced  by  the  geometrical  theory  of  diffraction  is  a 
cylindrical  wave  emanating  from  the  edge  of  the  wedge  and  is  de¬ 
scribed  by  a  diffraction  coefficient  obtained  from  the  solution  to  the 
perfectly  conducting  wedge  problem. 
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For  the  two-dimensional  case  illustrated  in  Fig.  10,  the  far- 
field  diffraction  for  a  unit  cylindrical  wave  incident  upon  a  perfectly 
conducting  wedge  is  given  as  a  sum  of  geometrical  optics  (v'  )  and 
diffracted  (  Vg)  rays  as[  5] 

(18)  u(  4)  =  v*  (  x0,  £-<*)  -  vv  (  xQ,  £+a  )  +  vB(  x0,  n,  £-a) 

-  vB(  x0,  n,  i+a  )  , 

where  the  choice  of  signs  is  determined  by  the  appropriate  boundary 
condition  (  -  for  u=0  and  +  for  the  normal  derivative  of  u=0  on  the 


walls) 

and  the  variables  are  defined 

are  given  by[  5] 

exp  [  jkr  cos 

(  19) 

v'"  (  r,  <t>)  -  • 

0,  otherwise 

and 

(20) 

Vg(  r,  n,  ()>) 

=  1T’3  e  ?  (- 

N=  0,1,2,... 


(» sin  ») 


2  [  cos(4>/2)[ 


T T  4> 

COS  ~  ~  COS  — 

n  n 


-i  t2 


ikr  cos  6  (  -1*.  j 

eJ  J  6  dT  l 

[  kr(  1+  cos4>)  ]2 


+  [  higher  order  terms] 


where  r  and  4>  are  dummy  variables.  For  large  values  of 
[  kr(  1  +  cos  41)  ]  >  vg  can  be  expressed  as 

7 T 

1 


-  J"?  -  jkr 

e  e 


(21) 


vB(  r,  n,4>)  - 


^2tt  ki 


IT  s*n  n 


4> 


cos  —  -  cos  — 
n  n 


X  [1  +  (  Mkr)"1  +  ...] 
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Fig.  1 0-- Three- dimensional  diffraction  by  a 
wedge  of  finite  edge  length. 


Using  the  first  term  in  the  asymptotic  expansion  given  in  Eq.  (21) 
write  the  diffraction  coefficient  for  plane  wave  incidence  as 

e-jf  e-Jkr 


(22) 


D(  r,  n,  41)  = 


\J  2tt  kr 

The  total  diffracted  far  zone  field  for  plane  wave  incidence  is  then 


TT  I 

i  sin 

7T 

cos  n  - 

c°s  n 

we 


(  23) 


uD(  r,  n,  4)=[  D(  r ,  n,  <j>")  -  D(  r,n,c|>+  )  ]  uincident 
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A  detailed  discussion  of  the  diffraction  coefficient  is  given  in  Appendix 

IL 

The  diffraction  of  a  spherical  wave  by  an  edge  has  been  studied 
by  Oberhettinger[  13]  and  also  by  Nomura[  14]  .  Oberhettinger  has 
obtained  the  diffraction  by  an  infinite  wedge  for  a  point  source  utilizing 
a  Green's  function  technique.  Nomura  has  obtained  the  solution  for  dif¬ 
fraction  by  an  infinite  wedge  for  a  dipole  source.  However  the  two- 
dimensional  diffraction  solution  can  also  be  applied  to  wedge  diffrac¬ 
tion  of  a  spherically  incident  field  through  the  use  of  ray  optic  tech¬ 
niques.  That  is,  the  magnitude  of  the  incident  ray  is  determined 
using  the  spatial  attenuation  factor  for  a  point  source  rather  than  that 
of  a  line  source,  and  the  diffracted  field  in  the  plane  of  the  point 
source  and  the  normal  to  the  wedge  is  then  obtained  using  Eq.  (  23). 

The  determination  of  the  diffracted  field  in  other  planes  requires  an 
extension  of  the  two-dimensional  solution  to  the  three-dimensional 
case.  Keller  has  shown  that  the  extension  of  Fermat's  principle  to 
edge  diffraction  in  three  dimensions  results  in  the  multiplication  of 
the  diffraction  coefficient  by  a  factor  of  l/sin  y  to  account  for  the 
dispersion  of  energy  in  the  cone  of  rays  illustrated  in  Fig.  9* 
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In  addition  the  source  of  rays  need  not  be  isotropic,  but  may 
have  a  field  pattern  in  the  general  case.  Thus  in  the  application  of 
Eq.  (  23)  to  diffraction  by  a  source  with  a  nonuniform  pattern,  the 
geometrical  optics  terms  are  multiplied  by  the  pattern  and  the  dif¬ 
fraction  terms  are  multiplied  by  the  value  of  the  pattern  in  the  di¬ 
rection  of  the  edge[  15]  .  Thus  the  far-field  form  of  the  diffraction 
of  a  source,  incident  at  an  angle  y  with  respect  to  the  edge  of  a 
perfectly  conducting  wedge,  is  given  by 


(  24) 


E(£) 


p<e) 

sin  y 


{  v'  (  XQ,  t  -a)  } 


+ 


p(ZTr-e) 

sin  y 


{v*  (xQ,  £+a)  } 


+  -P--^0>  (vB(  xQ,  n,  i-a)  -  vB(  x0,  n,  i+a)  }  , 

sin  y 

where  P(  £)  is  the  pattern  of  the  line  source  and  “  tt  +  a  is  the 
direction  of  the  edge.  The  phase  reference  for  Eq.  (  24)  is  the  edge 
of  the  wedge.  The  pattern  of  the  reflected  term  is  P(  2tt-£).  The 
factor  1  /sin  y  expresses  the  effect  of  conical  diffraction  and  xQ  is 
the  minimum  distance  from  the  source  to  the  edge. 

For  two-dimensional  diffraction  the  foregoing  technique  yields 
diffracted  rays  which  are  parallel,  having  been  derived  for  an  infinite 
edge.  This  diffraction  solution  is  not  valid  for  a  wedge  having  an  edge 
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of  finite  length*  However  it  will  be  assumed  that  the  diffraction  pat¬ 
tern  in  the  plane  containing  the  finite  edge  is  the  same  as  the  radi¬ 
ation  pattern  of  a  line  source  of  finite  length.  The  line  source  has 
an  excitation  determined  by  the  diffraction  by  an  infinite  wedge.  This 
assumption  takes  into  account  the  distribution  of  energy  in  three 
dimensions  of  the  family  of  parallel  diffracted  rays  due  to  an  edge 
of  finite  length.  The  effects  of  the  finite  edge  upon  the  radiation  pat¬ 
tern  in  the  plane  containing  the  edge  are  then  calculated  by  consider¬ 
ing  a  travelling  wave  on  the  edge  having  a  phase  velocity  of  (  kQ  cos  y) . 
Ignoring  end  effects,  the  radiation  pattern  in  the  plane  containing  the 
edge  is  that  of  the  travelling  wave  antenna  of  length  i  and  phase 
velocity  (  kQ  cos  y) .  This  pattern  can  be  evaluated  using  the  radiation 
integral.  Thus  the  total  diffraction  for  a  finite  wedge  shown  in  Fig.  10 
can  be  approximated  as 

k  i  sin  X 

(25)  E(0,4)  =  E(£)  -f™  - 

2  Xq 

where  E(  4)  is  the  excitation  obtained  using  the  diffraction 

coefficient  ,  and 

k0i  (cos  0  -  cos  y) 

x°= - 7 -  ' 

This  technique  has  been  applied  successfully  in  the  computation  of  the 
radiation  pattern  of  a  rectangular  waveguide  using  diffraction  tech- 
niques[  16]  .  The  use  of  the  diffraction  coefficient  in  combination 
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with  the  radiation  integrals  allows  the  determination  of  the  scattered 
field  in  three  dimensions  for  a  finite  wedge. 

B.  Diffraction  by  a  Pair 
of  Wedges 

The  process  of  diffraction  by  a  pair  of  interconnected  wedges  is 
shown  in  Fig.  11.  It  is  assumed  that  an  incident  plane  wave  strikes 
the  edge  of  each  wedge  of  included  angles  ax  and  a2.  The  incident 


Fig.  ll--Single  diffraction  by  a  pair  of  connected  wedges. 
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rays  are  diffracted  giving  rise  to  the  singly  diffracted  rays  and 

Uq2.  These  rays  may  be  evaluated  using  the  techniques  previously 
described.  In  order  to  calculate  higher-order  effects  an  extension 
of  the  geometrical  theory  of  diffraction  is  applied.  Namely,  dif¬ 
fraction  for  cylindrical  wave  incidence  is  applied  to  determine  the 
interaction  between  the  wedges. 

Since  the  incident  waves  are  plane  waves  the  plane  wave  dif¬ 
fraction  coefficient  is  used  to  calculate  the  singly  diffracted  rays. 
Taking  wedge  1  as  a  phase  reference  we  may  write 


(  26) 


UDi<0>  =  ^  sin  £ 
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7T  01-  0 

cos  - -  -  cos - 

nl 
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+  TT  Z7T  -  01  -  0 

COS - COS  - 

ni 


and 


(27) 


1  7 T 

UD ,(  °)  =  —  sin  — 
UZ  n2  n2 


COS 


7T 
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0-01 
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1 


TT  0  +  01 

COS -  -  COS  - - — 

n2  n2 


COS  0 


where  (  Z-rii  )  ir  =  orj  ,  (  Z-n2)  tt  =  az,  S  =  i  cos  01  and  the  reference 
for  the  angles  4>-,4>+  is  the  common  surface  of  the  two  wedges. 


32 


The  factor  e  ^  77  ^  2ukr  is  omitted  as  only  angular 

variations  are  of  interest.  When  the  angle  of  scattering  (  0)  is 
equal  to  the  angle  of  incidence  (  0)  as  in  the  case  of  backscatter, 
Eqs.  (  26)  and  (  27)  become. 
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(29) 


TT  -  1  .  TT 
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Figure  12  illustrates  the  shadow  boundary  and  the  reflected  waves 
which  exist  for  plane  wave  incidence  upon  the  pair  of  wedges.  In 
the  directions  of  the  reflected  fields  and  in  the  direction  of  the 
shadow  boundary  care  must  be  taken  in  the  evaluation  of  the  total 
field.  In  these  directions  a  geometrical  optics  field  exists,  and 
the  combination  of  the  terms  representing  the  geometrical  optics 
field  and  the  singly  diffracted  field  to  obtain  the  total  field  is 
specified  in  Appendix  II.  The  behavior  of  the  singly  diffracted 
fields  for  the  case  of  backscatter  when  the  incident  angle  0 
approaches  tt /Z  is  of  interest  as  it  is  seen  that  the  expressions 
for  the  singly  diffracted  fields  approach  infinity  as  0  approaches 
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Fig.  12  - -Geometrical  optics  rays  for  a  pair 
of  connected  wedges 

TT  /2.  If  the  singly  diffracted  fields  are  written  in  the  form  of  real 
and  imaginary  parts,  the  combination  of  Eqs.  (  28)  and  (  29)  re¬ 
sults  in 


(  30) 


Real  (  Ujj  +  Up>2)  =  cos  (  ki  cos  0) 
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(31) 


Imaginary  (Up^  +  Uj-^)  =  s*n  ^  ^  COS 
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When  the  real  and  imaginary  parts  are  placed  over  least  common 
denominators  and  the  limits  taken  through  application  of 
LfHospitalf  s  rule  it  is  found  that  the  real  part  diverges  while  for 
£  =  0  -  7 r  /2  small  the  imaginary  part  becomes 

(  32)  Imaginary  (Up  +  Up^)  =  sin  (  ki  sin  C) 
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Taking  the  limit  as  t,  “**  small,  Eq.  (  32)  becomes 


(33)  Imag  (  UDl  +  U^)  =  + 
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1  sin(  ki  sin  £,) 
2nx  sin  (  £/ni  ) 

1  sin (  ki  sin  £,) 
2ni  sin(  C/n2) 


+  ki  for  £  =  0 
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Keller[  1]  has  examined  this  problem  of  determining  the  scattered 
field  in  a  geometrical  optics  region.  He  presents  a  solution[  l]  using 
the  "Cross  Section  Theorem[  l]  n  or  "Forward  Scattering  Theorem!  15]  11 
that  states  that  the  forward  scattering  cross  section  Q  is 
(34)  Q  =  4 it  /k  Imag(  f(  0)  ) 

where  f(  0)  is  the  scattered  field  amplitude  in  the  forward  (  i.  e.  ,  geo¬ 
metrical  optics)  direction.  Noting  that  the  geometrical  optics  scat¬ 
tered  fields  exists  only  in  the  forward  and  back  directions,  then 
Im  f(  0)  is  the  geometrical  optics  result.  Using  this  theorem  Keller[  l] 
demonstrates  that  the  cross  section  obtained  using  single  diffraction 
may  be  obtained  in  this  way.  He  further  states  that  if  the  geometrical 
optics  terms  are  included  in  the  evaluation  of  the  scattered  field  the 
singularities  of  the  singly  diffracted  field  on  the  shadow  boundaries 
are  cancelled  resulting  in  a  finite  expression  for  the  scattered  field. 

The  behavior  of  the  total  field  on  the  shadow  boundary  is  discussed  in 
Appendix  IL  The  proper  combination  of  the  diffracted  and  geometrical 
optics  fields  results  in  continuity  of  the  field  across  the  shadow  bounda¬ 
ry. 

The  cross  section  theorem  expressed  by  Eq.  (  34)  yields  the 
geometrical  optics  cross  section  if  the  imaginary  part  of  the  singly 
diffracted  field  is  known.  This  imaginary  part  of  the  field  is  obtained 

TT 

e  "J4 

using  the  result  of  Eq.  (  33)  and  re-introducing  the  factor  — — =  from 

V  2TTk 

the  diffraction  coefficient.  The  backscattered  field  in  the  specular 
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direction  by  a  body  of  revolution  having  the  profile  shown  in  Fig,  11, 
may  be  obtained  using  the  above  results  and  the  spatial  attenuation 
factor.  The  geometrical  optics  backscattered  field  in  the  specular 
direction  may  thus  be  written  as 


(  35) 


F(  r) 


where  F(  r)  is  the  spatial  attenuation  factor.  An  approximate  solu¬ 
tion  for  the  spatial  attenuation  factor  for  a  conical  generator  is 
presented  in  Appendix  III.  In  the  region  close  to  the  specular 
direction  the  sin  (  x)  /x  pattern  behavior  indicated  by  Eq.  (  33)  may 
be  used. 

C.  Multiply  Diffracted  Rays 

Another  class  of  rays  which  are  of  interest  are  the  multiply- 
diffracted  rays.  An  example  of  multiple  diffraction  may  be  found 
in  the  analysis  of  parallel-plate  waveguides  by  Ryan  and  Rudduck[  1 8] 
and  in  coupling  between  parallel- plate  waveguides  by  Dybdal,  Rudduck 
and  Tsai[  1 9 ]  .  Figure  13  illustrates  the  doubly-diffracted  rays 
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Fig.  13- -The  multiply  diffracted  rays  for  a 
pair  of  connected  wedges. 
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UD2Di  (  6)  and  Uj-)  xj-)2(  0)  .  These  rays  arise  when  a  singly  diffracted 
ray  is  incident  upon  a  wedge  and  is  itself  diffracted.  In  the  case  de¬ 
picted  in  Fig.  13  the  doubly  diffracted  rays  may  be  expressed  by  as¬ 
suming  that  the  singly  diffracted  rays  are  cylindrical  rays  emanating 
from  the  edge  of  the  wedge.  Using  the  wedge  diffraction  coefficient 
for  cylindrical  fields  we  may  write 

(  36)  UD2Dl  (  0)  =  UDz(  0)  [  VB(  i  ,  nj  ,  O  +  vB(  i  ,  m  ,  <j>?  )  ] 

and 

(  37)  UDlD2^0)  =  UDi^  t  Vb(  1  '  n2»  ‘HF)  +  VB(  1  »  n2»  )  ] 

1  P  1  m  _ 

where  9l  =  9i  =  tt  -  0 

xP  —  i  nn  _  p 

<t>2  “  4>1  -  0 

For  spacings  between  the  wedges  on  the  order  of  two  wavelengths  or 
less  the  doubly  diffracted  rays  are  significant  as  has  been  demon¬ 
strated  by  Ryan  and  Rudduck[  16]  ,  and  Dybdal,  Rudduck  and  Tsai[  19]  . 
It  is  noted  that,  just  as  the  singly  diffracted  rays  correct  for  the 
discontinuities  in  the  geometrical  optics  field  which  arise  in  the 
directions  of  the  shadow  boundary  and  reflection  boundaries,  the 
doubly  diffracted  fields  correct  for  the  discontinuities  in  the  singly 
diffracted  fields.  Thus  it  is  seen  in  Fig.  13  that  the  singly  diffracted 
fields  Up)1(  0)  and  Up2(  0)  have  shadow  boundaries  at  0  =  tt  and  0 
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respectively.  This  process  continues  ad  infinitum,  that  is,  the 
triply-diffracted  rays  correct  discontinuities  in  the  doubly-diffracted 
rays  and  so  forth,  Yu  and  Rudduck[  19]  have  formulated  a  "Higher- 
Order  Diffraction"  technique  which  allows  the  effects  of  all  orders  of 
diffraction  to  be  calculated.  This  self-consistent  behavior  was  first 
described  by  Karp  and  Russek[  20]  .  This  technique  is  a  self-con¬ 
sistent  field  technique  which  generates  a  system  of  linear  equations. 
This  process  can  be  illustrated  using  Fig.  13,  where  it  is  assumed 
that  the  total  field  diffracted  at  each  wedge  can  be  expressed  as  the 
singly  diffracted  field  plus  the  diffraction  field  due  to  all  higher 
orders  of  diffraction  in  the  form 

(38)  UDl  (  6)  =  (  6)  U1  +  VB(  i  ,  m  ,  )  UD2(  0) 

and 

(39)  UD2(0)  =  D2(0)Ux  +  VB(i  ,n2,4»2)UDl(T7) 
where 

(40)  4>j  =  tt  -  0 
<j>2  =  0 

In  Eqs.  (  38)  and  (  39)  the  known  quantities  are  the  diffraction  coef¬ 
ficients  Di  (  0)  ,  D2(  0)  ,  VB(  i  ,  nj  ,  4>!  )  ,  VB(i,n2,4>).  Setting  0  equal 
to  tt  in  Eq.  \  38)  and  setting  0  equal  to  zero  in  Eq.  (  39)  results  in 
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(  41) 
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A  set  of  two  equations  in  the  two  unknowns  (  tt)  ,  and  U 0)  is 

determined,  thus  a  solution  is  obtained  for  the  total  diffraction  pat¬ 
tern 


(45)  Utot(0)  =UDl(0)  +UD2(6) 

This  form  of  solution  can  be  extended  to  the  interaction  of  any 
number  of  wedges.  In  most  cases  of  interest  it  is  not  necessary 
to  proceed  past  the  evaluation  of  the  doubly-diffracted  fields  in 
order  to  obtain  a  satisfactory  solution.  However  the  higher  order 
diffraction  solution  is  nearly  as  simple  as  the  evaluation  of  the 
doubly-diffracted  rays  and  is  thus  the  best  technique  in  general. 
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D.  The  Effects  of  Edge  Curvature 
on  Diffracted  Rays 

In  the  examples  of  diffraction  above  only  the  case  of  a  straight 
edge  was  discussed.  In  order  to  apply  the  wedge  diffraction  coef¬ 
ficient  to  more  general  cases  it  is  necessary  to  describe  the  behavior 
of  the  incident  and  diffracted  rays  for  a  curved  edge.  This  has  been 
done  by  DeVore  and  Kouyoumjian[  21  )  and  Ufimtzev[  22]  for  the 
circular  disk.  It  is  assumed  that  in  the  case  of  the  circular  disk 
shown  in  Fig.  14,  the  points  of  diffraction  which  contribute  to  the 
scattered  field  in  the  x-z  plane  are  located  at  cj)  -  0,tt.  It  is  further 
assumed  that  the  diffraction  coefficient  for  the  infinite  edge  given 
previously  can  be  assumed  to  be  valid  at  points  1  and  2.  This  as¬ 
sumption  is  consistent  with  the  property  of  diffraction  being  a  local 
phenomenon.  The  effects  of  curvature  on  the  behavior  of  reflected 
rays  has  been  described  previously.  This  analysis  can  be  extended 
to  diffracted  rays  using  the  extension  of  Fermat’s  principle.  Refer¬ 
ring  to  Fig.  15  we  have,  after  KouyoumjianJ  23] 


(  46)  c os  (3  =  i  •  t 

A  A 

-  sin  0'=  i  •  n 

A  A 
=  d  *  n 


0  <  (?  <  7T 


-  -  <  0*  < 
2  ~ 


3tt 

~T 


cos  6 


-  7T  <  6  <  7T 


42 


(b) 


Fig.  14- -Diffraction  by  a  circular  disk. 
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Fig.  15--Vector  relations  for  diffraction 
by  a  curved  surface. 

AAA 

where  t,  u,  b  are  the  tangent,  normal,  and  binormal  vectors  of  the 
edge  with  the  tangent  and  normal  vectors  in  the  plane  of  the  screen, 
and  where  i,  cl,  are  the  directions  of  incidence  and  diffraction.  The 
half  cone  angle  of  the  diffracted  rays  about  the  positive  tangent  is  (3. 
The  distance  p  between  the  edge  caustic  and  the  second  caustic  is 
given  by  Keller[  l]  as 

pI  =  -PQ  sin2p  /(  pQP  sin  |?  +  cos  \) 


(47) 
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Kouyoumjian[  2  3] 


(48) 


has  given  this 

a  •  (va> 

Pq  sin^  P 


relation  in  the  form 


where  pQ  is  the  radius  of  curvature  of  the  edge  at  Q. 

In  the  case  of  the  circular  disk  shown  in  Fig.  14,  the  caustic 
distances  pi  ,  p2,  can  be  obtained  from  Eq.  (48).  For  an  incident 
plane  wave  i  j  becomes  infinite,  and  for  backscatter  p  =  n/Z . 

Thus 

^  1  _  sin  0  +  sin  0^  1  „  sin  0  +  sin  0^ 

p2  a  p2  a 


The  negative  value  of  p2  means  that  the  caustic  is  between  the  edge 
and  the  field  point.  Thus  the  scattered  rays  pass  through  the  caustic 
and  a  phase  jump  of  tt/2  is  expected  to  occur  at  the  caustic.  This  is 
automatically  accounted  for  by  the  negative  value  of  p2. 

Equation  (48)  in  combination  with  ray  optics  thus  allow  the 
determination  of  the  effects  of  edge  curvature  upon  the  diffracted  rays. 
Assuming  that  the  edge  diffraction  coefficient  may  be  applied  at  each 
point  along  the  edge  the  diffracted  fields  of  a  curved  edge  may  be 


calculated. 
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E.  Correction  for  Axial  Caustics 

In  the  case  of  diffraction  by  a  circular  aperture  or  circular 
disk  for  an  axially  incident  plane  wave  each  point  on  the  edge  con¬ 
tributes  to  the  scattered  field  on  the  axis.  Thus  an  infinity  of  dif¬ 
fracted  rays  converge  on  the  axis  creating  a  caustic.  This  effect 
is  also  present  for  the  case  of  diffraction  by  a  slope  discontinuity 
on  a  body  of  revolution.  In  the  case  of  the  circular  aperture  the 
axis  is  also  the  direction  of  forward  scatter  and  for  the  case  of  the 
circular  disk,  the  axis  is  the  direction  of  specular  reflection.  The 
forward  scatter  theorem[  19]  allows  the  determination  of  the  fields 
in  the  geometrical  optics  regions  if  the  diffracted  field  can  be  ob¬ 
tained.  Keller[  l]  has  derived  a  ,fcaustic  correction  factor11  using 
an  asymptotic  solution  of  the  scalar  wave  equation.  By  comparison 
of  this  solution  with  the  divergent  result  calculated  using  the  dif¬ 
fraction  coefficient  and  ray  optics,  he  has  obtained  a  correction 
factor  which,  when  multiplied  by  the  diffraction  solution  yields  the 
correct  solution.  This  correction  factor  for  a  caustic  on  the  z-axis 
is 


(  50) 
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where  p  =  r  sin  0.  If  the  solution  is  desired  on  the  caustic  (  p  =  0)  , 
the  diffraction  solution  is  first  multiplied  by  the  correction  factor 
and  the  limit  as  p  tends  to  zero  is  taken. 

Another  procedure  for  obtaining  the  field  on  the  axis  in  the 
geometrical  optics  region  of  the  circular  aperture  or  disk  is  to  use 
the  equivalent  physical  optics  fields  over  the  aperture  (or  on  the 
disk)  and  evaluate  the  radiation  integral  for  these  fields.  The  case 
of  interest  for  backscatter  is  that  of  backscatter  by  a  circular  disk 
of  radius  a,  which  results  in  the  well-known  formula 


ka2  T  2  J i  (  2ka  sin  0) 


Es  =~- 


(  51) 


cos  0 


(  2ka  sin  0) 


The  case  of  axial  backscatter  by  a  slope  discontinuity  on  a 
body  of  revolution  can  be  treated  using  a  combination  of  diffraction 
theory  and  the  radiation  integral.  The  diffraction  coefficients  are 
used  to  specify  the  diffraction  at  each  point  on  the  "ring11  and  the 
radiation  integral  is  used  to  sum  these  contributions.  In  order  to 
obtain  the  individual  contributions  at  each  point  the  incident  wave  is 
decomposed  into  tangential  electric  and  magnetic  fields  in  the  x-y 
plane  at  the  edge  of  the  ring  shown  in  Fig.  1  6a.  Taking  the  plane 
of  incidence  to  be  the  x-z  plane  we  have  for  the  }  incidence  plane 


-  cos  0  sin  <j>!  1 


jka  sin  0  cos  cj)1  i 
r  e  JEl 


(  52) 


+  COS  <\>' 


47 


Fig.  l6--Coordinates  for  the  ring  source 


48 


/ 


(53) 


f  \ 

f 

-  COS  4> 1 

jka  sin  9  cos  6'  i 

E  plane 

-  cos  0  sin  4)' 

^  * 

-  e  H 

r 

H  plane 
v  > 

Applying  the  diffraction  coefficient  to  these  tangential  components  of 
the  incident  field  we  obtain  for  the  scattered  fields  at  the  edge  of  the 
ring 


(54) 


>- 


-  cos  6  sin  (V 
cos  4>  1 


k  [  vB(  r,  n,  ijj  ')  -  vB(  r,  n,  i|j  +  )  ]  E1 
jka  sin  0  cos 


•  e 


(  55) 


-  cos  <j> 1 
-cos  0  sin  cj)1 


■  [  Vg(  r,  n,  v|j  ')  +  Vg(  r,  n,  i|>+)  ]  H1 


Jka  sin  0  cos  <\>' 


where 


4^ 


A 

Vi 

A 

n 


s 


-  77  +  2s#  cos-1  (  -Vi  •  n) 

=  0 

=  incident  unit  vector 
=  unit  normal  to  the  surface 

+  1  if  -Vi  X  n  is  -y  directed 

A  A  A 

-1  if  -Vi  X  n  is  +y  directed 


The  scattered  fields  obtained  using  the  diffraction  coefficients  can  be 
related  to  equivalent  electric  and  magnetic  currents  at  the  edge  of  the 
ring  using  the  asymptotic  form  of  the  diffraction  coefficients  and  the 
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far  fields  of  electric  and  magnetic  current  filaments[  2 4]  .  For  z 
directed  filamentary  currents  we  have 


n 

-J* 


(  56) 


(  57) 


Ez  =  -z0kf 


Z\l  2-rrkp 


.7T 


Hz  =  -  Yck  lm  -JL^±_  .-J* 


Z^Z-nkp 


where  ZQ  and  Y0  are  the  impedance  and  admittance  of  free  space 
respectively.  Equating  Eqs.  (  56)  and  (  54)  ,  (  57)  and  (  55)  ,  and 


using  the  asymptotic  form  of  the  diffraction  coefficients  given  in 
Eq.  (  22)  results  in 

re  _  _2_  r  —  /  ,  -  v  —  /  « 

Zk 


(  58) 


r  =  -  ~  [G(n,40  -G(  n,  +)  ] 


''-cos  0  sin  1 


COS 


Jka  sin  0cos  <t>‘  ^i 


(  59) 


I1"  -  ~  [  G(  n,  ")  +G(  n,  +  )  ] 


-COS  4>' 


-cos  0  sin 
jka  sin  0  cos  i 


.JL 


where  D(  r ,  n,  i)j  ) 


“J4  -jkr 

2  e  J 

\l  2tt kr 


G(  n,  i|i ) 


These  equivalent  currents  on  the  edge  of  the  ring  apply  to  the  diffracted 
fields  only  and  do  not  include  geometrical  optics  terms. 
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Richmond[25]  has  presented  the  radiation  integrals  for  an  e- 
lectric  current  loop  of  radius  a  as 


(  60) 


2tr 

E0=  -  e  y  Ie(  <t>')  sin( 


jcopacos  0  -jkr. 


4ir  r 


0 


ejka  cos(  sin  6 


2ir 


(61)  E.-^e^of  Ie(<t>')  cos(  <^>-<t>') 

Y  4tjt~  J 


Jka  cos(<t>-4>1)  sin  6 


We  may  apply  the  Duality  Theorem  to  obtain  the  radiated  fields  in 
terms  of  the  magnetic  current  ( Im)  rather  than  in  terms  of  the  e- 
lectric  current  (  Ie)  ,  with  the  result 


2ir 


(62) 


H/ 


_  jco€  a  cos  6  -jk 


4nr, 


e-j^ro  y  Im(4>1)  sin(  <j>-<j>') 


jka  cos(  <J>-4>')  sin  6 


(63) 


=  -yoe4> 

TT  =  J^a 

H<t>  "  4irr0 


2tt 

e  ^r°  ^  Im(  <j>')  cos(4>-<t>') 
0 


gjka  cos(  )  sin  0 


d<)>' 


=  YOE0 
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In  the  case  of  backscatter  close  to  the  z-axis  0  =  01  with  0  small,  and 
we  take  4j+  =  'rr-2cos"1(z-n)  (  i.  e.  ,  0=0)  in  the  diffraction  coef¬ 
ficients.  Substituting  the  edge  currents  given  by  Eqs.  (  58)  and  (  59) 
into  the  radiation  integrals  we  determine  the  backscattered  fields. 
After  some  manipulation  we  obtain  for  the  E-plane, 


(64) 


E0=  -jaE1 


,•  „-jkr 


-  -cos20[G(n,iJj')  -G(n,iJj  +  )] 


Ji(u) 


+  [  G(  n,  iJj  )  +  G(  n,  ij;  +)  ] 
and  for  the  H -plane, 

(  65)  Hj  -  -jaE 


pTi(u) 


-J2(  u) 


■  jkr 


+  cos20[  G(  n,  ij;  )+G(n,  iJj  +  )]  ^  ^ 


u 


-  [  G(  n,  iJj  ')  -G(  n,  ij;  +  )  ] 


Ji  (  u) 


u 


J2(  u) 


where 


u  =  2ka  sin  0 


For  backscatter  on  the  z-axis  (  i.  e.  ,  on  the  Caustic)  we  have 


(  66) 


Eq  =  -  jaE'  G(  n,  ijj  ' ) 


+  v  e 


-jkr 


Thus,  in  the  case  of  backscatter  Eqs.  (  51)  and  (  66)  allow  the 
computation  of  the  axial  scattered  field  due  to  a  ’’ring11  slope 
discontinuity.  A  similar  integration  could  be  performed  for  non¬ 
circular  slope  discontinuities. 
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If  the  z-axis  is  a  geometrical  optics  region  we  must  apply  the 
cross  section  theorem[  l]  in  order  to  obtain  the  axial  scattered  field, 
and  thus  determine  the  equivalent  currents  on  the  ring.  In  order  to 
do  this  we  combine  the  integrals  of  Eqs.  (  60)  and  (  63)  to  obtain  for 
the  E- plane, 


(67) 
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2iTr 


Eg  =  -  Ja~  e  \  '|cos20  sin  <j)'  sin(4>-4>') 


1 
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IT  ,  tt  TT+211 

cos  —  —  1  cos  -cos 

n  n  n 
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IT  TT  TT  +  2  L 

COS  —  -1  cos - cos - — 
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^jka  sin  Geos  <j>!  ejka  sin  0  cos( 


where  -  sin_1(  sin  Geos  <V) 


Expressing  the  diffracted  fields  in  terms  of  real  and  imaginary  parts 
and  applying  the  cross  section  theorem[  l]  results  in 


2tt 
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Taking  the  limit  as  0 -*■  0,  and  performing  the  integration  results  in 

(69)  e|  =  ^L2S±^ 

L  ro 

Which  agrees  with  the  physical  optics  on-axis  result  for  a  disk  given 
by  Eq.  (  51).  The  physical  optics  solution  is  accurate  within  3  dB 
for  a  /  X  >  0.  5[  23]  . 

The  scattered  field  predicted  for  a  ring,  as  given  by  Eq.  (  bb)  , 
may  be  in  error  if  higher-order  diffraction  terms  are  significant. 
However  these  terms  can  be  evaluated  and  their  contribution  computed 
in  the  form  of  Eq.  (  66) . 

A  particular  case  of  interest  is  where  several  “rings"  contri¬ 
bute  to  the  near-on-axis  field.  Such  a  case  is  illustrated  by  a  cylin¬ 
der,  or  by  a  conically  capped  cylinder.  For  the  case  of  the  cylinder 
the  specular  return  from  the  flat  end  would  be  dominant,  and  the 
contribution  of  the  rear  "ring"  can  be  neglected.  In  the  case  of  the 
conically  capped  cylinder,  the  contribution  from  the  cylinder  end 
can  be  significant  with  respect  to  the  contribution  from  the  "ring" 
formed  by  the  cone -cylinder  junction.  In  such  a  case  the  total  field 
is  the  sum  of  the  individual  contributions  of  the  "rings".  In 
evaluating  the  contribution  from  each  ring  the  integration  limits 
must  define  the  illuminated  region.  Thus  in  the  case  of  the  rear 

,  TT 

"ring"  of  a  cylinder  these  limits  would  be  -  7  .  For  "rings"  of 
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radius  a  separated  by  a  length  i  as  illustrated  by  Fig.  16  b  the  scat¬ 
tered  field  is  then 


(70) 


e-j2W 


cos  9 


E 


s 

0,  Z 


where  the  subscripts  1  and  Z  refer  to  the  front  and  rear  rings  re¬ 
spectively  and  E^  .  is  defined  by  Eq.  (  64)  ,  and  the  phase  reference 
0,  l 

is  at  the  front  "ring11  . 

If  the  on-axis  field  is  desired,  the  limit  of  Eq.  (  70)  as  0  — *•  0 
may  be  taken.  Equivalently,  the  incident  field  specified  by  dif¬ 
fraction  along  the  shadow  boundary  (  i*  e#  ,  2-  E*)  may  be  used  and 

Z 

the  integration  limits  extended  to  0-Ztt.  In  the  second  case  it  is 
necessary  to  evaluate  the  scattered  field  at  a  small  angular  distance 
from  the  axis  in  order  to  avoid  higher  order  diffraction  from  the 
front  nringn  , 

If  the  radii  of  the  rings  differ,  as  for  a  conical  frustrum,  the 
contribution  of  the  second  ring  must  be  evaluated  using  an  integration 
over  the  illuminated  portion  of  the  ring. 


CHAPTER  IV 

CREEPING  WAVE  ANALYSIS  OF 
BODIES  OF  REVOLUTION 


A.  The  Creeping  Wave  Concept 

The  concept  of  creeping  waves  was  introduced  by  Franz  and 
Deppermanf  26,  27]  for  the  interpretation  of  the  scalar  solution 
for  diffraction  by  a  circular  cylinder  or  a  sphere.  Senior  and 
Goodrich[28]  have  obtained  a  representation  similar  to  that  of 
Franz  and  Depperman  through  the  application  of  the  Watson  trans¬ 
formation  to  the  Mie  series  solution  for  the  sphere. 

Kouyoumjian[  29]  has  presented  a  creeping  wave  solution  for  the 
sphere  which  includes  all  higher-order  modes. 

A  general  illustration  of  the  creeping  wave  format  is  shown 
in  Fig.  17.  An  incident  plane  wave  is  diffracted  at  a  point  of 
tangency  (designated  A)  on  the  target.  A  portion  of  the  diffracted 
energy  is  trapped  at  this  point,  resulting  in  a  wave  which  propa¬ 
gates  on  the  surface  of  the  target,  shedding  energy  by  radiation 
as  it  progresses.  Finally,  this  wave  reradiates  at  B  in  the  scat¬ 
tering  direction  of  interest.  This  "creeping  wave"  can  thus  be 
described  by  diffraction  coefficients  at  the  points  of  diffraction  and 
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Fig.  17--General  concept  of  the  scattered  field 
due  to  creeping  waves. 
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reradiation,  by  an  attenuation  factor  to  account  for  radiation  losses, 
and  by  a  description  of  the  ray  path  geometry  on  the  target  traversed 
by  the  creeping  wave.  Thus  we  write,  a  general  form  for  the  creeping 
wave  scattered  field  as 


B 


ES  =  Da  Dr  G(  s)  e  ~  A 

CW  Ad 


(71) 


s 


where 


is  the  creeping-wave  scattered  field 


E 


cw 


is  the  diffraction  coefficient  at  A 


Dg  is  the  diffraction  coefficient  at  B 


y(  s)  is  the  creeping-wave  propagation  factor 


is  the  arc  length  along  the  path 


s 


G(  s)  is  the  convergence  factor  for  the  surface  rays. 


The  primary  task  is  the  determination  of  these  diffraction,  attenu¬ 
ation  and  ray  path  factors  for  a  general  body.  However  this  does 
not  appear  feasible.  Thus  it  is  necessary  to  evaluate  these  factors 
for  canonical  targets  whose  exact  solutions  are  available,  i.  e.  ,  the 
cylinder  and  sphere.  In  order  to  obtain  a  more  general  solution 
for  these  factors  it  is  necessary  to  utilize  experimental  data  to 
obtain  an  empirical  solution  for  more  general  targets  such  as  the 
prolate  spheroid.  It  is  the  purpose  of  this  chapter  to  examine 


both  the  exact  and  empirical  solutions  for  the  sphere  in  order  to 
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formulate  a  solution  for  the  more  general  target.  It  will  be  demon¬ 
strated  that  a  rather  simple  approach  yields  accurate  results  for  a 
large  class  of  targets. 

The  form  of  the  simplified  creeping  wave  analysis  desired  is 
such  that  the  convergence  factor  expressed  in  Eq.  (  71)  may  be  taken 
equal  to  unity,  i.  e.  ,  no  convergence  or  divergence  of  rays  need  be 
considered.  In  addition,  it  is  desired  to  evaluate  the  creeping  wave 
contribution  due  to  all  creeping  waves  by  evaluation  of  Eq.  (71)  for 
a  single  equivalent  creeping  wave.  The  simplified  analysis  thus 
uses  a  single  equivalent  non-divergent  (or  non-conver gent) 
creeping  wave  to  construct  the  approximate  solution  for  the  scat¬ 
tered  field.  In  order  to  demonstrate  the  validity  of  the  simplified 
creeping  wave  solution  it  is  necessary  to  identify  the  major  creeping 
wave  contributors  in  existing  creeping  wave  solutions.  Consequently 
the  surface  fields  of  a  target,  for  which  the  exact  surface  fields  are 
known,  has  been  examined.  In  particular,  the  solution  of  interest  is 


that  for  the  sphere. 
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B.  The  Creeping  Wave  Solution 
for  Scattering  by  a  Sphere 

The  analysis  of  the  scattering  by  a  sphere  has  been  one  of  the 
most  important  problems  in  electromagnetic  theory.  The  classical 
solution  of  Mie  using  separation  of  variables  and  series  techniques 
is  the  basic  starting  point  for  casting  the  exact  solution  into  the 
asymptotic  form  on  which  the  creeping  wave  format  is  based.  The 
analysis  of  Senior  and  Goodrich[  2  7]  is  an  illustration  of  this  adaption* 
A  recent  paper  by  Hong[  30]  gives  a  derivation  which  is  more  closely 
related  to  the  geometrical  properties  of  the  ray  paths,  deriving  the 
attenuation  and  diffraction  coefficients  based  upon  the  properties  of 
the  ray  paths  from  the  exact  solution.  Hong’s  analysis  yields  the 
higher-order  correction  terms  for  the  cylinder  and  sphere  but  is  not 
applicable  to  more  general  bodies.  Hong  restricts  the  ratio  of  radii 
of  curvatures  in  the  propagation  and  orthogonal  directions  to  be  less 
than  or  equal  to  unity.  This  restriction  limits  his  analysis  to  the 
cylinder  and  sphere.  Keller  and  Levy[  31,  32]  have  also  treated  this 
problem*  Kinberf  33]  has  developed  a  general  technique  using  a 
M  semi-geodetic"  ray  coordinate  system  which  applies  to  bodies  of 
revolution  and  has  used  this  technique  to  develop  a  creeping  wave 
analysis  of  the  spheref  34]  .  Moreland,  Peters,  and  Kilcoyne[  34] 
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have  adopted  a  different  approach  which  is  empirical  in  nature  and 
utilizes  a  simplified  ray  path  geometry.  These  solutions  and  their 
adaption  to  more  general  targets  will  be  discussed  briefly. 

The  analyses  of  Kouyoumjian[  29]  and  Hong[30  ]  are  based  upon 
the  concept  that  the  incident  energy  is  diffracted  at  the  shadow  bounda¬ 
ry,  and  a  portion  of  this  energy  becomes  attached  to  the  surface  of 
the  sphere.  This  trapped  (  or  creeping)  wave  traverses  the  sphere 
on  a  geodesic  (  a  great  circle  route)  specified  by  the  direction  of  the 
incident  ray  shedding  energy  tangentially  as  it  progresses,  finally 
contributing  to  the  backscattered  field  as  depicted  in  Fig.  18.  The 
great  circle  paths  are  assumed  to  intersect  the  back  point  of  the 
sphere,  forming  a  caustic  at  that  point.  A  consequence  of  this  ray 
path  picture  is  that  in  the  directions  of  axial  scatter  a  line  caustic 
is  formed.  In  the  analysis  of  Keller  this  requires  the  evaluation  of 
a  ncaustic  correction  factor”  for  axially  scattered  fields.  Also  this 
picture  does  not  account  for  the  magnitude  of  the  field  in  the  minor 
creeping  wave  plane  as  noted  by  Kazarinoff  and  Senior[36]  ,  In 
each  of  these  analyses  the  results  for  the  echo  area  of  the  sphere 
are  in  good  agreement  with  the  exact  solution.  The  diffraction 
and  attenuation  coefficients  derived  by  Keller  and  Levy[  31,  32]  are 
carried  out  to  first  order  while  the  forms  given  by  Hong[  30]  and 
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Fig.  18- -Creeping  waves  on  the  sphere. 

Voltme r[  37]  incorporate  higher-order  correction  terms.  A  com¬ 
parison  of  the  results  of  Keller  and  Levy[  31,  32]  and  Voltmer[  37] 
is  given  in  Table  I  which  has  been  prepared  by  Voltmer[  3  7]  . 

An  approach  to  the  general  formulation  of  the  creeping  wave 
paths  on  bodies  of  revolution  has  been  developed  by  Kinber[  33]  . 
Kinber  has  shown  that  the  wave  equation  may  be  expanded  in  a  set 
of  "ray  coordinates"  in  which  the  solution  can  be  written  in  a  form 
whose  magnitude  is  dependent  upon  the  cross  section  of  the  ray  tube 


DIFFRACTION  AND  ATTENUATION  COEFFICIENTS 
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and  whose  phase  is  dependent  only  upon  the  propagation  constant  and 
the  path  length  traversed.  In  particular  he  seeks  a  solution  of  the 
form 

(72)  u  =  eikeW(£-C,  r,-0 

where  £,  r|  are  the  ray  coordinates.  In  the  case  of  the  sphere  the 
ray  coordinate  system  obtained  by  Kinber  [34]  is  shown  in  Fig.  19. 


Fig.  19--The  ray  coordinate  system  developed  by  Kinber. 
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Here  it  is  seen  that  the  ray  path  does  not  intersect  the  back  point  but 
is  instead  tangent  to  a  cone  whose  cone  angle  can  be  written  as 


(  73)  pQ  =  sin"1  (  m/ka) 


for  a  source  field  having  an  e  azimuth  dependence.  The  cone 


angle  is  dependent  upon  the  mode  order  (  m) .  The  analysis  of  Kinber 
for  an  arbitrary  incident  field  requires  a  summation  of  modes.  As 
the  ray  coordinates  used  for  each  mode  differ  the  task  of  constructing 
the  surface  field  at  a  given  point  on  the  sphere  is  formidable.  For 
the  case  of  plane  wave  incidence  of  the  form 


(74) 


the  analysis  for  the  scattered  field  is  simplified  resulting  in  a  solution 


for  the  scattered  field  of  the  form 


where  p,m  are  electric  or  magnetic  dipole  moments 

respectively 


h  is  the  distance  to  the  observation  point 


Y  =  kh/M,  Y!  =  kh'/M 
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This  solution  is  similar  to  that  of  Kouyoumjian[  29]  and 
Hong[  30 ]  ,  the  primary  difference  being  in  the  ray  path  geometry 
used.  A  complete  discussion  of  this  solution  is  given  by  Kinber[  34]  . 

To  date  the  above  solutions  have  not  been  applied  to  the  task  of 
determining  the  surface  fields  on  the  sphere.  Kazirinoff  and 
Senior[  36]  have  attempted  to  apply  creeping  wave  theory  to  calcu¬ 
late  the  major  and  minor  axis  surface  fields.  While  the  agreement 
in  the  major  axis  is  good,  the  minor  axis  fields  have  not  been 
satis  fa  ctorialy  calculated. 

The  empirical  approach  adopted  by  Moreland,  Peters,  and 
Kilcoyne[  35]  utilizes  the  known  E-plane  field  to  determine 
appropriate  attenuation  and  diffraction  coefficients.  This  approach 
is  of  interest  in  that  it  lends  itself  to  extension  to  more  general 
targets.  It  proceeds  by  assuming  that  the  surface  field  shown  in 
Fig.  18  in  the  shadow  region  of  the  sphere  can  be  represented  as 
the  sum  of  two  creeping  waves  as 
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(  76) 


H*  =H0 


(  jk0  +  as)  a(  tt/2  -  0) 


+  e 


-( jkQ  +  «s)  a(  tt/2  +  0) 


and  the  radial  electric  field  may  be  written  as 


(  77) 


:r  =  ed  j^_(jko  +  as>  a(*/2  -  0) 
.'(  jkQ  +  as )  a(  tt/2  +  0)J 


-  e 


Now  HQ  and  E0  are  chosen  to  match  the  exact  solution  as  closely  as 
possible  for  some  value  of  kQa.  A  value  of  kQa  =  10  was  chosen. 
Then,  assuming  as  to  be  modified  by  a  constant  from  the  value  of  that 
given  by  Keller  and  Levy[  32]  for  the  cylinder,  one  obtains 


(  78) 


2  JL  4?* 

ag  =  0.  84  a  3  X” 3  e^ 


in  order  to  fit  the  curves  of  Fig.  20. 

It  is  seen  that  fields  of  the  form  given  in  Eqs.  (  76)  and  (  77) 
are  a  reasonable  approximation  to  the  exact  fields.  The  back- 
scattered  fields  due  to  the  creeping  wave  may  be  written  as 


(  79) 


ES  =  2E: D2  e 
cw  i  s 


-j2k0a  e~(  jkQ  +  as)  ira  e  jk°R 

R 


where  R  is  the  observation  distance  from  the  center  of  the  sphere. 
The  magnitude  and  phase  of  the  creeping  wave  obtained  using  Eq. 

(  79)  can  be  compared  to  the  magnitude  and  phase  of  the  exact 
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ANGLE  B 


Fig.  20- -Normalized  fields  on  the  surface  of  a 
conducting  sphere. 
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Fig.  2l--Absolute  relative  magnitude  of  backscattered 
fields  due  to  creeping  waves. 
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creeping-wave  fields  obtained  by  subtracting  the  specular  point  contri¬ 
bution  of  the  physical  optics  term  from  the  exact  scattered  fields;  i.  e.  , 


(80) 


exact 

'CW 


=  E 


exact 


(  koa) 


(' 


1 

j2k0a 


) 


a 

2 


where  the  factor  e  /r  has  been  suppressed.  The  value  of  D2  is 

adjusted  by  a  constant  to  agree  with  the  exact  backscattered  fields  for 

2  i 

a  radius  of  1.  OX.  The  form  of  Dg  is  also  modified  by  a  factor  of  X  in 
order  that  Eq.  (  79)  be  dimensionally  correct,  resulting  in 


1  1  _;tt/12 

(  81)  D2  =  0.  27  p3  X3  e  J 


The  magnitudes  of  the  exact  and  approximate  creeping-wave  com¬ 
ponents  as  a  function  of  radius  are  compared  in  Fig.  21  and  their 
phases  are  compared  in  Fig.  22.  The  agreement  is  good  for  a  wide 
range  of  kQa  values. 

Moreland,  Peters,  and  Kilcoyne[  35]  have  also  applied  this 
form  of  solution  to  the  problem  of  bistatic  scattering  in  the  E-plane 
with  excellent  results. 

The  analysis  of  Moreland  et.  al.  [35]  is  a  simplified  creeping 
wave  analysis.  It  uses  a  single  non-convergent  (and  non- divergent) 
ray  path  together  with  approximate  diffraction  and  attenuation  co¬ 
efficients  to  construct  an  approximate  solution  for  the  scattered  field 
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due  to  creeping  waves  propagating  in  each  direction  along  the  ray  path. 
The  path  chosen  is  the  path  traversed  by  the  "major"  ray  (i.  e.  the  path 
corresponding  to  the  E-plane  of  the  sphere).  This  analysis  suggests 
that  an  approximate  picture  of  scattering  by  a  sphere  can  be  con¬ 
structed  by  neglecting  the  creeping  waves  which  have  a  radial 
magnetic  field  (i.  e.  the  "minor"  creeping  waves)  and  by  considering 
only  the  creeping  waves  which  have  a  radial  electric  field  (i.  e. 
the  "major"  creeping  waves). 

The  analyses  outlined  above  do  not  allow  the  determination 
of  the  fields  in  the  minor  axis  because  of  the  ray  geometry 
which  has  been  assumed.  The  analysis  of  Kinber  [34]  may  allow 
such  a  calculation  but  this  has  not  yet  been  performed. 

A  new  approach  to  the  determination  of  the  ray  paths  and  the 
surface  fields  on  the  sphere  is  considered  here.  It  will  be  assumed 
that  all  the  modes  may  be  lumped  into  a  single  traveling  wave  or 
that  only  a  single  mode  is  dominant.  In  order  to  determine  the  path 
followed  by  such  a  creeping  wave  the  results  of  the  exact  Mie  series 
solution  are  used  to  compute  the  real  part  of  the  Poynting  vector 


at  the  shadow  boundary,  i.  e.  , 


(DEGREES) 
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Fig.  22--Phase  of  backscatter ed  fields  due  to 
creeping  waves. 


(  82) 
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4> 


It  is  assumed  that  a  creeping  wave  originating  at  a  point  on  the 
shadow  boundary  (  tt  /  Z ,  4> )  propagates  in  the  direction  of  Real  (S)  at  that 
point  and  thereafter  follows  a  great  circle  route  around  the  sphere. 
The  result  of  such  an  assumption  is  shown  in  Fig.  23  where  it  is 
seen  that  the  ray  paths  do  not  intersect  the  back  point  of  the  sphere, 
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Fig.  23--Ray  paths  on  the  sphere  determined 
using  the  Poynting  Vector. 
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and  indeed  miss  the  back  point  by  varying  amounts.  In  addition  to 
the  ray  paths  the  creeping  wave  fields  have  been  calculated  along 
the  ray  path  in  the  form 

-(  jk_  +  Qs)i 

(83)  Hcw  =  (H0cos4>)W2*  e  J°  *F j[  i  ) 

4> 

where  the  factor  cos  4>  accounts  for  the  tangential  component  of  the 
incident  magnetic  field  at  the  shadow  boundary  and  where  the  exact 
as  is  used  and  i  is  measured  along  the  great  circle  specified  by  the 
Poynting  vector.  The  calculations  were  performed  using  a  computer 
program.  For  the  case  of  the  ka  =  10  sphere  the  intersections  of  the 
ray  paths  in  the  <f>  s  0°,  90°  planes  have  been  plotted  in  Figs.  24  and 
25.  In  addition  the  angles  at  which  the  ray  paths  cross  the  axes 
are  plotted  in  Figs.  26  and  27.  As  seen  from  Fig.  23  the  H-plane 
magnetic  field  can  be  written  as 

(84)  H0  =  Hq  cos  4>  F(  i  )  e"^k°  +  1  sin  6 

where  6  is  the  angle  between  the  ray  and  the  H-plane  and  where 
the  attenuation  constant  as  is  the  attenuation  constant  for  the 
"hard"  sphere  given  in  Table  L  Using  the  intersection  and  angle 
data  from  Figs.  25  and  27  it  is  possible  to  determine  F(f  ),  l,  and  6 
for  each  ray  which  intersects  the  minor  axis,  and  thus  compute  the 
magnetic  field  in  the  H-plane.  The  results  of  such  a  calculation 
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Fig.  24--Intersection  of  the  ray  paths  and  the 
cj>  «  0°  plane  (i.  e.  ,  the  E- Plane). 
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Fig.  25- -Inter section  of  the  ray  paths  and  the 
4>  =  90°  plane  (i.  e.  ,  the  H-plane)* 
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Fig.  26--Angle  between  the  ray  path  and  the 
<f>  =  0°  plane. 
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boundary,  the  contribution  of  the  "minor"  creeping  wave  (  i.  e.  ,  the 
rays  corresponding  to  the  "soft"  sphere)  may  be  important,  even 
though  the  attenuation  of  this  wave  is  large,  as  can  be  seen  from 
Table  L  Using  the  attenuation  and  diffraction  coefficients  of  Table 
I,  and  assuming  no  convergence  or  divergence  of  the  minor  creep¬ 
ing  wave  about  the  <j>  =  90°  plane  the  minor  creeping  wave  may  be 
calculated.  The  results  of  this  calculation  are  also  shown  in  Fig. 

28  .  Close  to  the  shadow  boundary  the  minor  creeping  wave  field 
is  seen  to  be  a  reasonable  approximation  to  the  exact  field.  Thus 
both  the  major  and  minor  creeping  wave  contributions  are  needed 
to  predict  the  minor  axis  field.  A  more  accurate  solution  for  the 
effect  of  the  major  creeping  wave  near  the  shadow  boundary  is 
needed  in  order  to  demonstrate  the  continuity  of  these  solutions  in 
the  4>  =90°  plane. 

In  the  calculation  of  the  major  axis  field  (<j>  =  0°  plane)  it 
suffices  to  know  the  intersection  of  the  closest  ray,  as  indicated 
in  Fig.  24.  This  is  because  the  rays  removed  from  the  major 
axis  have  little  effect  on  the  field  in  the  major  axis  due  to  polari¬ 
zation,  flux  tube  width,  and  crossing  angle  effects.  The  calculated 
major  axis  field  is  shown  in  Fig.  29.  It  is  seen  to  be  in  good  agree¬ 
ment  with  the  exact  solution.  The  calculated  field  is  normalized  to 
the  exact  field  at  0  =  ~  ,  <\>  =  0.  A  disadvantage  of  the  Poynting 
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Fig.  27--Angle  between  the  ray  path  and  the 
4>  =  90°  plane. 
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for  a  ka  =  10  sphere  are  presented  in  Fig.  28,  where  the  field  is 
normalized  to  the  exact  field  at  0  =  y  ,  <f>  =  0  .  It  is  seen  that  close 
to  the  back  point  the  calculated  field  is  in  good  agreement  with  the 
exact  solution.  The  agreement  deteriorates  as  one  approaches  the 
shadow  boundary.  This  is  expected  as  the  attenuation  coefficient 
for  the  "hard"  sphere  derived  by  Voltmerf  37]  is  most  accurate 
for  the  deep  shadow  region.  The  region  close  to  the  shadow  bounda¬ 
ry  has  been  examined  by  Foch[  38]  .  Also,  close  to  the  shadow 


Fig.  28--H-plane  surface  fields  on  the  sphere. 
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boundary,  the  contribution  of  the  "minor**  creeping  wave  {i.  e.  ,  the 
rays  corresponding  to  the  "soft1*  sphere)  may  be  important,  even 
though  the  attenuation  of  this  wave  is  large,  as  can  be  seen  from 
Table  I.  Using  the  attenuation  and  diffraction  coefficients  of  Table 
I,  and  assuming  no  convergence  or  divergence  of  the  minor  creep¬ 
ing  wave  about  the  <J>  as  90°  plane  the  field  of  the  minor  creeping  wave 
may  be  calculated.  The  results  of  this  calculation  are  also  shown 
in  Fig.  28.  Close  to  the  shadow  boundary  the  minor  creeping  wave 
field  is  seen  to  be  a  reasonable  approximation  to  the  exact  field. 

Thus  both  the  major  and  minor  creeping  wave  contributions  are 
needed  to  predict  the  H-plane  field.  A  more  accurate  solution  for 
the  effect  of  the  major  creeping  wave  near  the  shadow  boundary  is 
needed  in  order  to  demonstrate  the  continuity  of  these  solutions  in 
the  <j>  a:  90°  plane. 

In  the  calculation  of  the  major  axis  field  (<f>  =  0°  plane}*  it 
suffices  to  know  the  intersection  of  the  closest  ray,  as  indicated 
in  Fig.  24.  This  is  because  the  rays  removed  from  the  major 
axis  have  little  effect  on  the  field  in  the  major  axis  due  to  polari¬ 
zation,  flux  tube  width,  and  crossing  angle  effects.  The  calculated 
E-plane  field  is  shown  in  Fig.  29.  It  is  seen  to  be  in  good  agree¬ 
ment  with  the  exact  solution.  The  calculated  field  normalized  to 
the  exact  field  at  0=  ,  <)>  s  0.  A  disadvantage  of  the  Poynting 
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Fig.  29-  -  E-plane  surface  fields  on  the  sphere. 


vector  approach  to  determining  the  ray  paths  is  that  the  exact  solu¬ 
tion  at  the  shadow  boundary  must  be  known.  Thus  it  is  not  possible 
to  extend  this  technique  to  targets  for  which  the  exact  solution  is  not 
available.  Also  the  ray  paths  must  be  recalculated  for  each  sphere 
radius  of  interest.  These  calculations  have  been  performed  for 
values  of  kQa  =  4,  0  and  10,  0.  Over  this  range  of  kQa  values  the 


8? 


ray  paths  do  not  differ  substantially,  as  may  be  seen  from  Table  II 
where  some  major  and  minor  axis  intersections  are  compared 

It  is  noted  that  the  attenuation  coefficient  developed  by  Moreland, 
et  ah  [  35]  for  the  sphere  does  not  yield  the  correct  major  axis  fields 
when  the  effect  of  the  convergence  of  the  flux  tube  to  form  a  caustic 
is  considered.  This  is  reasonable  as  the  solution  of  Moreland,  et 
al.  [  35]  includes  the  effects  of  the  caustic  and  integration  of  the 
shadow  boundary  contributions  in  the  determination  of  the  empirical 
attenuation  and  diffraction  coefficients.  However  good  results  are 
obtained  for  the  backscatter ed  fields  using  their  solution.  The  fact 
that  their  approximate  solution  gives  reasonable  results  is  not 
surprising  if  one  considers  the  ray  path  picture  of  Kinber[  34]  and 
that  of  the  Poynting  vector  method.  It  is  seen  that  in  the  E -plane 
the  major  axis  creeping  wave  is  the  dominant  contributor,  and  the 
rays  which  skew  about  the  back  point  are  negligible  as  regards  back¬ 
scatter  since  there  are  no  tangent  lines  from  the  skewed  creeping 
wave  paths  in  the  direction  of  the  source.  These  creeping  waves 
would  contribute  to  the  bistatic  radar  cross  section  in  planes 
removed  from  the  E-plane.  Thus,  because  of  the  ray  path  geometry, 
a  simple  traveling  wave  picture  suffices  to  describe  the  E-plane. 


MAJOR  AND  MINOR  AXIS  INTERSECTIONS 
FOR  A  SPHERE  (  ka  =  4.  0,  10.  0) 
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The  Poynting  vector  analysis  suggests  that  a  ’’simplified  ray 
path  geometry”  based  upon  the  component  of  the  electric  field  normal 
to  the  surface  may  be  used  for  more  general  targets.  In  such  a 
geometry  only  the  electric  field  component  normal  to  the  surface 
(  i.  e.  ,  the  ’’major  creeping  wave”)  is  considered.  If  the  family  of 
rays  is  considered,  thus  determining  the  spatial  attenuation  factor 
F(  i  )  the  exact  diffraction  and  attenuation  coefficients  of  Keller  and 
Levy[  27,  28]  or  Voltmer[  32]  are  used  to  compute  the  creeping  wave 
fields.  If  only  the  ray  corresponding  to  the  point  at  which  the  incident 
E-vector  is  normal  to  the  surface  is  used  and  a  single  nondivergent 
(  or  nonconvergent)  path  is  used,  the  approximate  diffraction  and  at¬ 
tenuation  coefficients  of  Moreland,  et  aL  are  used. 

A  simplified  ray  path  geometry,  using  the  approximate  forms  of 
the  attenuation  and  diffraction  coefficients  has  the  advantage  of  being 
easy  to  apply  to  the  case  of  parallel  polarization  backscatter  for 
complicated  targets  as  the  geodesic  path  is  simply  the  profile  of  the 
target.  In  the  case  of  perpendicular  polarization  the  geodesic  path 
is  more  difficult  to  determine.  A  numerical  method  for  finding  a 
geodesic  path  is  presented  in  Appendix  IV.  Several  examples  of  a 
simplified  analysis  will  now  be  considered  for  the  prolate  spheroid, 
the  ogive,  and  the  spherically  capped  ogive. 
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C.  The  Prolate  Spheroid 

The  prolate  spheroid  represents  a  body  for  which  an  exact 
closed  form  solution  is  not  available  except  on  the  axis  of  rotation. 

The  creeping  wave  solution  for  the  prolate  spheroid  proposed  here 
is  similar  to  that  of  the  sphere.  That  is,  it  is  assumed  that  the  inci¬ 
dent  plane  wave  is  diffracted  at  the  shadow  boundary  causing  the  exci¬ 
tation  of  creeping  waves  which  propagate  around  the  spheroid  and 
contribute  to  the  scattered  field.  The  scattering  due  to  the  illumin¬ 
ated  portions  of  the  spheroid  will  be  computed  using  geometrical 
optics. 

In  order  to  find  the  ray  path  geometries  for  the  prolate  spheroid 
for  an  arbitrary  angle  of  incidence,  the  geodesic  corresponding  to  the 
point  of  attachment  and  the  tangent  direction  at  that  point  must  be 
calculated.  As  seen  from  the  sphere  analyses  of  Kinber[  34]  and  by 
the  Poynting  vector  method  the  determination  of  the  tangent  direction 
of  the  creeping  wave  at  the  shadow  boundary  is  not  an  easy  task  in 
general. 

The  prolate  spheroid  has  previously  been  studied  by  Moffatt[39] 
using  time  domain  analysis.  Moffatt's  approach  is  similar  to  the 
creeping-wave  approach  in  that  he  approximates  the  impulse  response 
through  the  use  of  a  ’’physical  optics”  contribution  in  the  form  of  an 
impulse  at  time  t  =  0,  which  decays  and  joins  into  a  creeping  wave 
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return  at  time  t  =  TQ  (where  T0  corresponds  to  the  time  required  for 
the  creeping  wave  to  propagate  around  the  spheroid).  The  form  of 
the  impulse  response  postulated  by  Moffatt  is  shown  in  Fig.  30. 


Fig.  30--Impulse  response  of  the  prolate  spheroid. 


Moffatt1  s  determination  of  the  transit  time  T0  of  the  creeping 
wave  was  carried  out  by  finding  the  path  length  of  the  ellipse  defined 
by  the  intersection  of  the  plane  defined  by  the  incidence  vector  and 
the  incident  E  vector  with  the  prolate  spheroid.  The  initial  point  of 
intersection  is  the  point  at  which  the  surface  normal  and  the  E-vector 
are  collinear.  This  form  of  analysis,  which  resulted  in  excellent 
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agreement  suggests  that  the  simplified  ray  path  geometry  used  by 
Moreland,  Peters,  and  Kilcoyne[  35]  for  the  sphere  can  be  applied 
to  the  prolate  spheroid. 

Following  this  line  of  reasoning  we  write,  for  the  case  of 
parallel  polarization  for  the  prolate  spheroid  that 
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In  the  case  of  the  prolate  spheroid  the  creeping  wave  paths  are  el¬ 
liptical,  thus  the  total  attenuation  must  be  expressed  as  an  integral 
which  is  dependent  upon  the  radius  of  curvature  along  the  path.  Also 
the  radius  of  curvature  at  the  points  of  attachment  and  reradiation  (A,  B) 
must  be  computed  in  order  to  determine  the  diffraction  coefficient. 

It  was  found  that  the  product  of  the  diffraction  coefficients  at 
the  points  of  attachment  and  reradiation  could  be  related  to  the  dif¬ 
fraction  coefficient  of  the  sphere  with  the  sphere  radius  replaced  by 
\l  pg  .  The  square  of  the  diffraction  coefficient  is  thus 

(86)  D2s  =  0.  27  (  pApB)  X3  e  J 


where  and  are  the  radii  of  curvature  in  the  propagation  direction 
at  points  A  and  B  respectively.  The  specific  solution  for  the  prolate 
spheroid  has  been  presented  by  Ryan[  40]  and  Peters  and  Ryan[4l]  in 
previous  papers.  This  case  is  automatically  formulated  by  the 
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computer  program  for  the  second  degree  surface  of  revolution.  This 
program  compute  the  geometrical  optics  contribution,  the  points  of 
attachment  and  reradiation  and  performs  the  integration  over  the  path 
for  the  case  of  parallel  polarization. 

The  investigation  of  the  prolate  spheroid  carried  out  previously 
indicated  that  the  magnitude  of  the  attenuation  factor  required  for  ac¬ 
curate  results  in  the  case  of  the  prolate  spheroid  differs  from  that  of 
the  sphere.  Computations  were  performed  using  different  magnitudes 
and  the  results  are  shown  in  Fig.  31.  It  was  concluded  that  the  form 
of  the  approximate  attenuation  coefficient  appropriate  to  the  2:1 
prolate  spheroid  for  parallel  polarization  was 

(87)  a  =  0.  55  p"^  \"3 

where  p  is  the  radius  of  curvature  in  the  propagation  direction. 

Previous  work  by  Peters[  42]  had  been  devoted  to  the  determin¬ 
ation  of  the  effects  of  the  ratio  of  the  orthogonal  radii  of  curvature  on 
the  attenuation  and  diffraction  coefficients.  Using  the  case  of  the 
ogive  he  was  able  to  determine  that  the  diffraction  coefficient  was  the 
same  as  given  in  Eq.  (86)  but  that  the  attenuation  coefficient  was 
modified  to 

cxog  =  0.  20  p“*  \"3  eiz 


(88) 
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Fig.  3  l{a)--Backscatter  patterns  of  the  prolate  spheroid  for  a 
range  of  magnitudes  of  the  attenuation  coefficient. 
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(b) 


Fig.  3  l(b,  c)--Backscatter  patterns  of  the  prolate  spheroid  for  a 
range  of  magnitudes  of  the  attenuation  coefficient. 
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It  is  seen  that  the  attenuation  coefficient  for  the  2;  1  prolate 
spheroid  examined  lies  between  the  values  for  the  sphere  and  ogive. 
Figure  29  illustrates  the  magnitudes  of  the  approximate  attenuation 
factor  for  the  single  path  simplified  ray  geometry  as  a  function  of 
the  orthogonal  radii  of  curvature  for  bodies  which  have  been  examined 
to  date.  It  is  seen  that  the  points  may  be  fitted  by  the  smooth  curve 
shown  in  Fig.  32  of  the  form 


(89) 


a  =  ^1 


-0.  84(  pp/pn) 


Pp 


3  X"3  ej  6 


1.  48  e  "  ■'  hP/hO'  +  o.  20 

where  pp  and  pQ  are  the  radii  of  curvature  in  the  propagation  and 
orthogonal  directions  respectively. 


Fig.  32- -The  variation  of  the  magnitude  of  the  approximate 
attenuation  coefficients  as  a  function  of  the  two 
orthogonal  radii  of  curvature. 
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This  form  has  been  used  to  compute  the  scattered  fields  of  the  prolate 
spheroid  for  parallel  polarization  and  is  in  agreement  with  the  best 
fit  obtained  in  Fig.  31  using  Eq.  (87).  In  addition  the  echo  area  as 
a  function  of  wavelength  was  calculated  for  parallel  and  perpendicular 
polarization  for  0  =  90°  .  These  results  are  shown  in  Figs.  33  and  34 
and  are  seen  to  be  in  excellent  agreement  with  the  experimental  data 
of  Moffatt.  [  39  ] 


2M  PROLATE  SPHEROID 


Fig.  33--Echo  area  of  the  prolate  spheroid  as  a  function  of 
kQa  for  perpendicular  polarization  and  0  =  90°. 
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K  (SEMI-MINOR  AXIS) 


Fig.  34--Echo  area  of  the  prolate  spheroid  as  a  function  of 
kQa  for  parallel  polarization  and  0  =  90°. 
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The  form  of  the  attenuation  coefficient  given  by  Eq.  (  89)  yields 
good  results  when  used  in  conjunction  with  a  simplified  ray  path  ge¬ 
ometry  which  utilizes  a  single  ray  path  defined  by  the  "major11  creep¬ 
ing  wave.  The  results  of  Moffatt[  39]  indicate  that  this  form  of  the 
ray  path  geometry  is  also  applicable  to  the  case  of  perpendicular 
polarization.  Due  to  the  difficulties  of  computing  the  ray  paths  for 
the  general  body  for  perpendicular  polarization,  computations  have 
not  yet  been  made  for  this  case.  This  will  be  a  future  goal  in  the 
extension  of  the  creeping  wave  analysis. 

D.  The  Ogive 

The  ogive  has  been  treated  by  Peters[  43]  using  a  traveling 
wave  antenna  approach.  In  the  near-nose-on  region  this  solution 
has  been  the  only  satisfactory  solution  to  date.  The  form  of  this 
traveling  wave  solution  suggests  that  a  creeping  wave  solution  for 
the  ogive  can  be  found  providing  that  the  effects  of  the  tip  in  the 
shadow  region  can  be  determined.  Since  no  diffraction  coefficient 
for  a  cone  tip  is  available  an  analytical  expression  for  the  reflection 
and  transmission  coefficients  at  the  ogive  tip  is  not  available.  How¬ 
ever,  experimental  measurements  have  been  performed  for  a  range 
of  tip  angles  to  determine  reflection  coefficients.  Using  these 
measurements  Peters  [42]  has  determined  the  ogive  attenuation 


factor  previously  cited. 
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Figure  35  shows  the  ogive  and  the  coordinate  system.  As  in  the 
case  of  the  prolate  spheroid,  incident  energy  is  trapped  at  points  A 
and  B  where  the  incident  vector  is  normal  to  the  surface.  These 
creeping  waves  travel  along  the  ogive  until  they  encounter  the  tip. 


*- 

z 


Fig.  35--Coordinate  system  for  the  ogive. 


96 


At  the  tip,  part  of  the  creeping -wave  energy  is  diffracted  as 
creeping  waves  on  all  geodesics  containing  the  tip  and  the  remainder 
is  radiated  from  the  tip.  The  previous  treatment^  35]  assumed  that 
the  creeping-wave  fields  on  all  of  the  geodesics  have  equal  magni¬ 
tudes.  This  is  only  approximately  correct.  Small  differences  in  the 
magnitude  of  diffracted  rays  for  different  geodesics  will  yield  the  on- 
axis  creeping-wave  field  which  can  be  observed  in  the  measurements 
reported  by  Blore[  44]  .  Off-axis  the  creeping  wave  return  is  a 
result  of  the  creeping  waves  which  propagate  on  the  ogive  in  the  plane 
defined  by  the  incidence  direction  and  the  incident  E-field.  In  this 
plane  the  effect  of  the  tip  can  be  characterized  by  reflection  and 
transmission  coefficients  at  the  tip. 

Referring  to  Fig.  35,  the  backscatte red  fields  attributed  to  the 
creeping-wave  fields  may  be  written  as 

ES  =  E1  •  D2  {Rt[e"^k°+  a°g^  2Ll  e'2koS 

CW  t 

+  e'(  jko  +  aog)  ■  2L2  e^2k°S] 

-jkR 

+  2Tt[e’(  Jko+  aog)  (  Li  +  L2)  1  )  —7= — 
tl  6  J  }  \[2it  R 


(90) 
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whe  r  e 


D2  is  the  diffraction  coefficient  given  in  Eq.  (8  6) 

Rt  is  the  tip  reflection  coefficient 

is  the  tip  transmission  coefficient 
aQg  is  the  ogive  attenuation  coefficient 
Lj  ,  L2  are  shown  in  Fig.  35,  and 

S  refers  the  phase  to  the  center  of  the  ogive. 


For  the  spherical  ogive  Rj-  =  -T^  =  F  approximately,  and  values  of 
F  have  been  obtained  experimentally.  Bistatic  radar  cross  section 
measurements  of  the  ogive  have  shown  that  the  energy  incident  on  the 
tip  is  scattered  by  the  tip  along  each  geodesic  on  the  ogive,  i.  e.  , 
these  tip  -  scatte r ed  fields  show  no  <j>  dependence.  This  agrees  with 
the  results  of  the  equivalent  antenna  approach,  for  which  no  <f>  vari¬ 
ation  of  the  field  was  assumed.  Thus  a  factor  of  1  /\TZtt  is  introduced 
to  account  for  the  lack  of  <|>  dependence  of  the  tip-scattered  fields. 

The  radar  cross  section  of  the  ogive  calculated  using  Eq.  (90)  with 
F  =  0.  54  is  compared  with  the  experimental  results  and  with  the  echo  area 
calculated  by  the  equivalent  antenna  approach  in  Fig.  36.  In  this 
figure  the  equivalent  antenna  solution  has  been  corrected  for  radi¬ 
ation  loss  using  the  ogive  attenuation  factor  as  described  by 
Peters[  42]  . 
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Fig.  36- -Backscattered  field  patterns  for  the  ogive. 


Moreland,  et  al.  have  applied  creeping -wave  theory  to  analysis 
of  a  spherically  capped  ogive  as  shown  in  Fig.  37.  This  body  repre¬ 
sents  another  target  having  a  discontinuity  in  the  shadow  zone.  They 
have  expressed  the  scattered  field  as[  35] 


(91) 


E_  =  E.D2  e"2v°gr0°  {  T  [  e‘j2koi  +  ZyogrQ  +  e~ZyogrQ  ] 


+  2T  T 

os  so 


[g-Oko*  +Vsrs29l  ♦)»l|}  1  e-jkR 

n/2ttR 
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where 


V 


°g 


®o» 


D2 

Vs 


=  jko  +  aog  • 

-  radius  of  curvature  of  the  ogive 

-  radius  of  the  spherical  cap, 

=  angles  shown  in  Fig,  37, 

=  attenuation  coefficient  of  ogive, 

=  diffraction  coefficient,  given  in  Eq.  (8  6) 

=  Jko  +  as  ’  and 

=  attenuation  coefficient  of  the  sphere. 
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The  reflection  coefficient,  T,  and  the  transmibsion  coefficients,  Tos 
and  Tso,  arise  because  of  the  change  of  surface  impedance  at  the 
junction  of  the  ogive  and  sphere.  The  surface  impedance  of  the  ogive 
is  taken  to  be  ZQ,  the  impedance  of  free  space,  and  the  surface  im¬ 
pedance  of  the  sphere  is  taken  to  be  the  ratio  of  E0  to  H0  for  the  fields 
of  the  major  creeping  waves  postulated  in  Eqs.  (  76)  and  (  77)  for  the 
sphere.  Referring  to  Fig.  20  we  see  that  this  ratio  is  dependent 
upon  the  sphere  radius  and  can  be  determined  from  the  exact  fields 
of  the  sphere  for  a  given  radius.  Thus  F  and  TQS  can  be  written  as 


(92) 

and 

(93) 


r  - 


os 


z  s  Zp 
Zs  +  Zo 


2ZS 

Zs  +  zo 


It  is  necessary  to  make  T  the  negative  of  Eq.  (92)  in  order  to  obtain 
agreement  with  the  shape  of  the  experimental  curves.  This  sign  ap¬ 
pears  since  the  discontinuity  acts  a**  a  line  scatterer  or  a  caustic. 

The  ray  passing  through  it  undergoes  a  phase  shift  of  tt/2  radians  and 
the  transmitted  ray  passes  through  two  such  discontinuities.  In  ad¬ 
dition,  the  factor  1  / \^2tt  proposed  by  Moreland,  et  al.  is  retained 
to  account  for  the  rotational  symmetry  of  the  caustic  which  causes 


diffracted  rays  to  propagate  along  the  spherical  cap  in  directions 
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other  than  that  of  the  plane  of  incidence,  thus  decreasing  the  backscat- 
tered  field.  Combining  all  the  terms  the  total  expression  for  the  back- 
scattered  field  is 


(94) 


Es  =  EiD2e"Z7ogr0°  {rle'j^o1  + 2^ogr 0 +  e'2^ogr0] 


+  2  TosT,n[  e'^  ko£  +  ^srs2el  +  j2lT)  ]  }  __I -  e"jkR 

\/2ttR  R 


The  results  of  the  computation  for  the  spherically  capped  ogive  are 
presented  in  Fig.  38  and  compared  with  experimental  measurement. 
It  is  seen  that  in  the  near  nose-on  region  the  results  are  in  good 
agreement  with  measurement.  This  form  of  creeping -wave  analysis, 
coupled  with  the  traveling  wave  picture  and  physical  optics,  allows 
prediction  of  the  backscattered  field  of  the  ogive  over  a  wide  range 
of  incidence  directions  and  kQa  values.  It  should  be  noted  that  the 
impedance  of  the  creeping  wave  is  obtained  by  a  best  fit  of  the  fields 
at  the  surface  of  the  sphere  in  the  entire  shadow  region.  This  would 
include  the  fields  close  to  the  point  at  the  rear.  It  has  been  noted 
that  this  is  a  region  of  poor  fit,  regardless  of  the  approximation 
used,  because  there  are  fields  present  here  which  are  not  included 
in  the  simple  creeping-wave  picture.  This  creeping  wave  picture, 
as  noted  previously,  does  not  include  a  caustic  at  the  point  at  the 
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Fig.  38- -  Backscattered  field  patterns  of  the 
spherically  capped  ogive. 
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rear.  Further  improvement  in  these  results  would  be  anticipated  if 
the  fields  were  more  closely  matched  by  two  creeping  waves  at  a 
distance  removed  from  the  rear  point  and  if  multiple  interactions  are 
inc  luded. 

It  is  seen  that  the  effects  of  discontinuities  in  the  shadow  zone 
must  be  carefully  accounted  for.  Future  work  in  this  area  devoted 
to  the  determination  of  the  reflection  and  transmission  coefficients 
of  tips  and  wedge  discontinuities  would  provide  a  more  accurate 
solution  for  the  effects  of  shadow  zone  discontinuities. 

In  summary  the  attenuation  and  diffraction  coefficients  cor¬ 
responding  to  the  simplified  single  ray  geometry  are  presented  in 


Table  IIL 


ATTENUATION  AND  DIFFRACTION  COEFFICIENTS  FOR  THE 
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is  the  radius  of  curvature  in  the  direction  orthogonal  to  propagation. 
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Et  The  Extension  of  Creeping- 
Wave  Analysis  to  Disks 

The  simplified  ray  path  geometry  suggests  the  examination  of 
radar  targets  where  only  one  possible  creeping-wave  path  exists. 
Such  a  class  of  targets  are  disks,  for  incident  plane  waves  in  and 
polaiized  in  the  plane  of  the  disk.  A  new  solution  for  this  class  of 
targets  may  be  obtained  by  assuming  that  a  creeping-wave  field  can 
exist  on  the  edge  of  the  disk  in  such  a  case  and  that  the  functional 
forms  of  the  diffraction  and  attenuation  coefficients  which  have  been 
us lC  for  volumetric  shapes  apply.  Thus  we  write  for  one  com¬ 
ponent  of  tne  cr  ^eping-wave  field 


(  5) 


whu  e 


rcw  _  ^i  ^2  -<2ji  -j  -jkR/R 

E  =  E  e  a  e  ‘  e  J 

ECW  is  the  scattered  creeping-wave  field 

E1  is  the  incident  field 

is  the  complex  attenuation  coefficient 

f  is  the  path  length  traversed  on  the  edge 

D2  is  the  product  of  the  diffraction  coefficients 
d 


at  the  points  of  attachment  and  reradiation. 
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For  elongated  targets  such  as  ogives  the  radius  of  curvature  in 
the  direction  of  propagation  (  p)  of  the  creeping  wave  is  much  greater 
than  the  orthogonal  radius  of  curvature  (  pQ)  for  the  case  of  near¬ 
nose-on  backscatter  where  the  creeping  wave  contribution  is  domi¬ 
nant.  If  the  circular  disk  is  considered  to  be  a  limiting  case  of  an 
oblate  spheroid  as  the  ratio  of  minor  to  major  axes  becomes  small, 
it  is  seen  that  the  ratio  (  p/pQ)  is  large  as  in  the  case  of  the  ogive. 
Thus  it  is  reasonable  to  apply  the  ogive  attenuation  coefficient  to  the 
disk.  Thus  for  the  disk  we  take 


-jtt/6  iL  1 
a  =  0.  20  e~J  p"3  \~3 


(96) 


where  p  is  the  radius  of  curvature  in  the  direction  of  propagation. 


In  order  to  obtain  the  diffraction  coefficient  it  is  assumed  that  the 
functional  form  is 


2  1 


(97) 


DZd  =  A  p3  \3 


We  now  write  the  scattered  field  of  the  ogival  disk  as 


(98) 


Ecw  =  y  D2  E1  (  e"2(+jko  +  a)il  +  e"2* 
d 

-2  e-^  j^o  +  Q  )  (  -^l  +  ^2)  ) 


-2(  +Jko+  «)  *2 


where 


y  is  the  measured  voltage  reflection  coefficient 
of  the  tip  (+0,  73  for  the  ogive  treated)  and 


f  l  and  are  shown  in  Fig.  39. 


Now,  adjusting  the  constant  A  to  agree  with  data  as  shown  in  Fig.  39 
we  find  that  A  =  0.  1  approximately.  We  now  apply  the  above  attenu¬ 
ation  and  diffraction  coefficients  to  the  cases  of  circular  and  elliptical 
disks.  However,  we  first  need  to  obtain  an  expression  for  the  dif¬ 
fracted  field  at  the  specular  point  for  these  targets. 

For  a  disk  the  diffracted  field  at  the  specular  point  can  be 
written  using  Eq.  (  6)  as 


(99) 


Esp  =  EXDW  's/d  /(  d  +  R) 


where  Esp  is  the  specularly  diffracted  field 


D 


w 


is  the  wedge  diffraction  coefficient 


d  is  the  distance  from  the  specular  point  to  the 
caustic  of  the  diffracted  rays  in  the  plane  of 
the  disk  and 


R  is  measured  from  the  specular  point. 


In  the  case  of  the  circular  disk  shown  in  Fig.  40  we  have 


d  =  a/Z.  For  a  disk  the  wedge  diffraction  coefficient  becomes 


(  100) 


D. 


w 


-e 


-j(kR  +  ir/4) 
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THICKNESS  •  0.0625*  V**'4" 

(a) 


Fig.  39--Backscattered  field  pattern  of  the  ogival  disk. 
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Fig.  40--Scattering  by  a  circular  disk  edge-on. 


Combining  Eqs.  (99)  and  (  100)  and  taking  the  phase  reference  at 
the  specular  point  we  may  write  for  the  circular  disk 

(101)  Esp  =  -  I  E1  \fa  /irk  e'j^  e"jkR/R 

A  similar  expression  can  be  derived  for  the  elliptical  disk.  For  the 
circular  and  elliptical  disks  we  now  write  the  total  scattered  field  as 


Etot  =  Ef*  +  ESP 

tot 


(  102) 
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The  radar  cross  section  for  the  circular  disk  and  results  for  a  2:1 
elliptical  disk  obtained  using  this  approach  are  presented  in  Figs. 

41  and  42.  In  some  cases  the  measured  data  is  presented  as  a 
line  representing  the  maximum  and  minimum  echo  as  read  from  the 
measured  pattern.  For  the  elliptical  disk  it  is  seen  that  good  results 
are  obtained  for  1.  5  <  ka  <  5,  while  the  results  for  the  circular 
disk  are  good  for  ka  >  1. 

Thus  empirical  attenuation  and  diffraction  coefficients  are 
obtained  which  yield  results  in  good  agreement  with  measurements 
for  ogival,  circular,  and  elliptical  disks.  The  agreement  between 
theory  and  experiment  demonstrates  that  the  concept  of  the  creeping 
wave  may  be  applied  to  disks,  for  waves  incident  in,  and  polarized 


in,  the  plane  of  the  disk. 


Ill 


Fig.  41--Radar  Cross  section  of  a  circular 
disk  edge- on. 
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a. 


Fig.  42{a)--Radar  cross  section  of  a  2:1  elliptical 

disk  edge-on 
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k  •  (  SEMI  -  MINOR  AXIS  ) 


b. 

Fig.  42(b)--Radar  cross  section  of  a  2;1  elliptical 

disk  edge-on. 
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K  •  (  SEMI  -  MINOR  AXIS  ) 


c. 


Fig.  42(c)--Radar  cross  section  of  a  2;  1  elliptical 

disk  edge-on. 
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d. 

Fig.  42(d)--Radar  cross  section  of  a  2;1  elliptical 

disk  edge-on. 
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CHAPTER  V 

COMPUTER  PROGRAM  RESULTS 


The  theory  discussed  in  the  preceeding  text  has  been  applied 
to  the  determination  of  the  radar  cross  section  of  targets  of  revo¬ 
lution  through  the  development  of  two  computer  programs.  These 
programs  are  discussed  in  Appendix  V.  The  programs  have  been 
tested  for  horizontal  polarization  backscatter  for  circular  cylinders, 
cones,  double  cones,  conically  capped  cylinders,  prolate  spheroids, 
a  prolate  spheroid- sphere  combination  and  a  prolate  sphe roid- oblate 
spheroid  combination.  The  targets  are  shown  in  Fig.  43,  and  the 
corresponding  patterns  are  shown  in  Figs.  44  to  49.  The  agree¬ 
ment  between  the  measured  and  calculated  echo  area  patterns  is  an 
indication  of  the  accuracy  of  the  various  techniques  used  in  the 
computer  program.  The  solution  corresponding  to  each  region  of 
the  pattern  is  identified  in  Figs.  44  to  49.  The  results  for  each  of 
the  targets  will  be  discussed  individually  in  the  following  paragraphs. 
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Fig.  43--Test  targets  for  the  geometrical  optics  and  creeping 
wave  sections  of  the  computer  program. 
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A.  Circular  Cylinder 

The  measured  and  calculated  patterns  for  the  circular  cylinder 
are  presented  in  Fig.  44.  As  indicated  in  the  figure,  the  solution  for 
the  axial  caustic,  the  single  diffraction  solution,  and  the  solution  for 
the  specular  direction  have  been  used  to  calculate  the  echo  area  pat¬ 
tern.  For  the  region  close  to  the  axial  caustic  (ie,,  0  close  to  0° 

or  180°)  the  agreement  is  excellent.  The  magnitude  of  the  echo  area 
on  axis  is  within  1  dB  of  measurement.  In  addition  the  location  of  the 
first  null  of  the  calculated  and  measured  patterns  is  the  same,  although 
the  depths  of  the  measured  and  calculated  nulls  do  not  agree.  For  the 
region  close  to  the  specular  direction  (  i.  e.  ,  0  close  to  90°)  the  agree¬ 
ment  is  excellent,  being  within  1  dB  at  the  specular  direction.  Again 
the  position  of  the  first  null  of  the  calculated  and  measured  patterns  agrees . 
The  depths  of  the  measured  and  calculated  nulls  are  also  in  reason* 
able  agreement,  having  an  error  of  approximately  5  dB.  The  re¬ 
mainder  of  the  pattern  is  calculated  using  the  single  diffraction  solu¬ 
tion.  This  region  shows  a  slight  null  shift  between  the  calculated  and 
measured  patterns.  The  magnitudes  of  the  calculated  and  measured 
peaks  are  within  3  dB.  The  accuracy  obtained  for  the  circular  cylin¬ 
der  is  within  +  1  dB  of  the  accuracy  obtained  by  Ufimtsev[  22]  and  is 
the  same  as  the  accuracy  obtained  by  Bechtel  and  Ross[45]  using 
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Fig*  44- -Monostatic  echo  area  of  the  circular  cylinder. 
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singly  diffracted  components.  These  results  demonstrate  that  the 
computer  solution  is  valid  for  a  cylinder. 

B.  Cone 

The  measured  and  calculated  patterns  for  the  cone  are  presented 
in  Fig.  45.  The  measurements  were  performed  by  Keys  and 
Primich[  46]  .  As  indicated  in  the  figure,  the  solution  for  the  axial 
caustic,  the  single  diffraction  solution,  and  the  solution  for  the 
specular  direction  have  been  used  to  calculate  the  echo  area  pattern. 

A  diffraction  solution  for  the  cone  has  also  been  presented  by 
Bechtel[  47]  which  did  not  include  an  axial  caustic  correction  and 
consequently  is  in  error  for  axial  incidence.  The  measured  data 
is  shown  as  a  range  of  values.  The  calculated  pattern  is  seen  to  be 
in  good  agreement  with  the  measured  values  over  the  entire  range 
of  incidence  angles.  However,  it  is  noted  that  the  accuracy  of  the 
measurements  is  not  good,  having  a  variation  of  approximately  5  dB. 
Thus  it  is  difficult  to  specify  an  absolute  accuracy  for  the  solution 
by  comparison  to  the  data.  However,  the  agreement  obtained 
indicates  that  the  computer  solution  is  valid  within  the  accuracy  of 


the  measurements. 
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Fig.  45- -Monostatic  echo  area  of  the  cone. 
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C.  Double  Cone 

The  measured  and  calculated  patterns  for  a  double  cone  are 
presented  in  Fig.  46.  As  shown  in  the  figure  the  solution  for  the 
axial  caustic,  the  single  diffraction  solution  and  the  solution  for  the 
specular  region  have  been  used  to  compute  the  echo  area  pattern. 

As  seen  from  the  figure  the  agreement  close  to  the  axial  caustics 
(  i.  e.  ,  0  close  to  0°  or  180°)  is  excellent.  On  the  caustic  the 
measured  and  calculated  echo  areas  are  within  \  dB,  and  the  lo¬ 
cation  of  the  first  nulls  is  in  error  by  only  2°  .  The  peaks  of  the 
specular  solution  are  within  2  dB  of  the  measured  peaks.  The  re¬ 
mainder  of  the  calculated  pattern  is  not  in  good  agreement  with  the 
measurement,  although  the  characteristics  of  the  pattern  are  indi¬ 
cated  by  the  calculations.  Comparison  of  the  measured  and  calcu¬ 
lated  patterns  shows  that  nulls  which  appear  in  the  calculated  pat¬ 
tern  are  filled  in  the  measured  pattern.  In  addition,  there  is  a 
peak  at  6  =  90°  which  is  not  predicted  by  the  calculations.  These 
results  indicate  that  higher -order  diffraction  terms  are  significant 
for  this  target,  at  aspects  removed  from  the  axis.  Also  as  this 
target  is  small  in  terms  of  wavelength  (a  =  0.  667X)  higher-order 
interactions  are  to  be  expected.  Improved  results  would  be  antici¬ 
pated  if  the  scattering  from  the  tips  of  the  cones,  and  the  inter¬ 
actions  between  the  joint  and  the  cone  tips  were  included  in  the 
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solution.  In  order  to  accomplish  this  task  the  diffraction  coefficient 
for  a  cone  tip  is  needed.  Such  a  diffraction  coefficient  for  arbitrary 
angles  of  incidence  and  diffraction  is  not  known.  Thus  further  theo¬ 
retical  study  is  required  before  a  more  accurate  solution  for  this 
target  can  be  obtained. 

D.  Conically  Capped  Cylinder 

The  conically  capped  cylinder  represents  a  more  complicated 
target,  which  approximates  a  practical  target.  The  measured  and 
calculated  patterns  are  shown  in  Fig.  47.  As  indicated  in  the  figure 
the  solution  for  the  axial  caustic,  the  single  diffraction  solution  and 
the  solution  for  the  specular  region  have  been  used  to  calculate  the 
echo  area  pattern.  The  regions  close  to  and  on  the  axis  (  i.  e.  ,  0 
close  to  0°  or  180°)  are  in  agreement  within  2  dB.  The  specular 
region  (0  close  to  90°)  is  in  agreement  within  1  dB.  The  specular 
region  normal  to  the  cone  ^  0  close  to  45°)  is  not  in  good  agreement, 
having  several  nulls  in  the  calculated  pattern  which  do  not  agree.  It 
is  noted  that  these  nulls  occur  at  the  angles  where  the  solution  shifts 
from  a  diffraction  solution  to  the  specular  solution  and  vice  versa. 

The  peak  of  the  specular  solution  is  within  3  dB  of  the  measured 
values  in  this  region.  It  is  noted  that  the  diffraction  solution  is  within 
3  dB  of  the  measurements  except  for  the  region  close  to  the  normal 


to  the  cone.  As  the  error  is  due  to  the  fact  that  the  calculated  nulls 
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capped  cylinder. 
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in  this  region  are  filled  in,  this  error  may  be  attributed  to  higher- 
order  diffractions  between  the  cone  tip  and  the  junction  between  the 
cone  and  cylinder,  and  to  scattering  by  the  cone  tip,  both  of  which 
have  been  neglected  in  this  solution.  An  additional  source  of  error 
may  be  attributed  to  the  target  used  for  the  measurement.  The  joint 
between  the  cone  and  cylinder  was  machined  with  a  radius  curve. 
Thus  the  wedge  diffraction  coefficient  may  not  accurately  describe 
the  actual  scattering  mechanism  at  the  joint.  In  spite  of  these  dis¬ 
agreements  the  over-all  character  of  the  scattering  pattern  is  pre¬ 
dicted  by  the  computer  program,  and  except  for  the  specular  region 
of  the  conical  section  the  average  error  is  less  than  3  dB. 

E.  Prolate  Spheroid 

The  solutions  used  for  this  smooth  target  are  the  geometrical 
optics  solution  for  the  specular  point,  and  the  creeping  wave  solu¬ 
tion.  As  a  test  of  the  general  creeping  wave  program  the  2:1  prolate 
spheroid  patterns  previously  computed  using  a  special  case  program 
which  were  presented  in  Fig.  31  were  re-computed  using  the  general 
program.  The  results  were  identical  to  the  results  presented  in 
Fig.  31  for  an  attenuation  coefficient  magnitude  of  0.  55.  Thus  the 
error  is  the  same  as  shown  in  Fig.  31,  being  less  than  3  dB  for  a 
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minor  axis  greater  than  \  wavelength*  These  results  indicate  that 
the  generalized  portions  of  the  computer  solution  were  operating 
correctly.  Next,  some  composite  shapes  were  examined. 

F.  Prolate  Spheroid -Sphere 
Combination 

The  prolate  spheroid- sphere  combination  represents  a  two 
section,  continuously  curved  target  with  no  discontinuity  in  the  first 
derivative  of  the  surface  at  the  join.  The  measured  and  calculated 
patterns  for  this  target  are  shown  in  Fig.  48.  As  in  the  case  of  the 
prolate  spheroid,  the  geometrical  optics  and  creeping  wave  solution 
were  used  for  all  aspect  angles  to  calculate  the  scattered  fields. 

The  agreement  between  measured  and  calculated  patterns  is  good, 
being  within  3  dB  for  aspects  more  than  30°  removed  from  the  axis. 
Within  approximately  30°  of  the  axis  of  revolution  the  accuracy  is 
not  as  good,  having  an  error  at  the  null  of  the  calculated  pattern  of 
5  dB.  Away  from  the  null  the  error  is  less  than  3  dB.  This  inac¬ 
curacy  close  to  the  axis  is  attributed  to  the  effects  of  the  discontinuity 
in  the  second  derivative  of  the  surface  at  the  join.  This  discontinuity 
results  in  a  reflection  at  the  joint  as  was  previously  discussed  in  the 
case  of  the  spherically  capped  ogive.  These  reflections  are  not  ac¬ 
counted  for  in  this  program.  Improved  results  would  be  expected  if 
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this  reflection  were  accounted  for.  An  additional  improvement  would 
be  expected  if  a  more  accurate  solution  for  the  scattered  fields  due 
to  the  illuminated  region  were  used. 

G.  Prolate  Spheroid-Oblate 
Spher ole  Combination 

The  measured  and  calculated  patterns  for  the  prolate  spheroid- 
oblate  spheroi combination  are  shown  in  Fig.  49.  This  target  is 
similar  to  the  prolate  spheroid- sphere  combinations  with  approxi¬ 
mately  the  same  results.  That  is,  accuracy  with  3  dB  for  aspects 
greater  than  40  removed  from  the  axis  of  revolution  with  errors 
as  large  as  8  dB  at  the  calculated  nulls  within  40°  of  the  axis.  The 
greater  discontinuity  in  the  second  derivative  of  the  surface  at  the 
join  for  this  case  results  in  a  larger  error  in  the  calculations  near 
the  axis  than  for  the  prolate  spheroid- sphere  combination.  Again 
the  addition  of  the  reflections  due  to  the  join  and  a  more  accurate 
solution  for  the  scattered  fields  due  to  the  illuminated  region  would 
be  expected  to  yield  improved  results. 

The  results  obtained  for  the  targets  discussed  above  show 
that  the  general  computer  program,  incorporating  first  order  dif¬ 
fraction  theory,  geometrical  optics,  and  creeping  wave  theory  is 
operable.  Improvements  in  this  program  have  been  suggested  by 
the  results,  and  include  the  use  of  multiple  diffraction,  improved 
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solutions  for  the  illuminated  region  and  creeping  wave  reflections. 

The  addition  of  multiple  diffraction  to  the  program  is  a  straight¬ 
forward  but  tedious  problem.  The  problems  of  an  improved  solu¬ 
tion  for  the  illuminated  region  and  of  the  reflections  of  the  creeping 
waves,  require  further  theoretical  study  before  they  can  be  included 
in  a  general  computer  solution. 

It  may  be  concluded  that  the  general  computer  program  yields 
results  which  are  useful  in  the  estimation  of  the  radar  cross  section 
of  composite  targets.  The  size  of  the  target  in  wavelengths  will  affect 
the  accuracy  of  the  solution.  The  restrictions  on  the  size  of  the 
target  for  a  class  of  targets  can  be  determined  by  comparison  of  the 
measured  and  calculated  values  for  a  special  case,  as  for  a  cylinder 
or  spheroid. 

In  general,  the  accuracy  of  the  computer  solution  will  improve 
as  the  size  of  the  target  in  wavelengths  is  increased.  Thus  if  a 
lower  limit  on  target  size  for  a  given  accuracy  is  established  for  a 
class  of  targets,  this  accuracy  will  be  maintained  for  larger  targets. 


CHAPTER  VI 
CONCLUSIONS 


The  geometrical  theory  of  diffraction  and  creeping  wave  solutions 
described  have  been  organized  in  computer  programs  to  calculate  the 
backscattered  fields  from  sectionally  continuous  bodies  of  revolution. 

The  results  obtained  indicate  that  these  techniques  are  applicable  to 
this  problem.  The  programs  have  been  used  to  investigate  targets 
having  either  straight  line  or  curved  smoothly- joined  profiles.  Targets 
having  combination  profiles  have  not  been  examined.  However,  the 
results  obtained  indicate  that  targets  having  such  combination  profiles 
can  be  treated  using  this  analysis,  and  that  an  extension  of  the  present 
computer  programs  to  treat  these  targets  is  feasible.  This  extension 
would  combine  the  existing  computer  programs  and  provide  the  ad¬ 
ditional  logic  required  to  treat  the  more  general  combination  target. 

The  computer  analysis  would  also  be  improved  by  including  the 
effects  of  interactions  (  i.  e.  ,  multiply  diffracted  waves  and  reflected 
creeping  waves)  .  The  theory  of  multiple  diffraction  has  been  described. 
The  reflection  effects  are  of  two  kinds:  reflection  from  tips,  and 
reflection  from  discontinuities  in  surface  derivatives  at  the  junction 
between  analytic  sections.  The  reflection  coefficients  for  tips  could 
be  obtained  experimentally,  and  the  results  fitted  by  an  empirical 
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curve.  The  effects  of  a  surface  derivative  discontinuity  can  also  be 
treated  experimentally,  although  physical  optics  is  known  to  give  a 
valid  result  for  this  contribution.  Thus  the  effects  of  reflections  could 
be  included  in  an  improved  computer  analysis. 

The  effects  of  fins,  ducts,  antennas,  or  other  special  scatterers 
may  be  treated  using  special  programming,  and  are  compatible  with 
the  existing  program.  Also  the  extension  of  this  analysis  to  concave 
targets  appears  feasible.  The  results  of  Hutchins  [ 5 0 ]  in  determining 
an  interior  wedge  diffraction  coefficient  could  be  applied  for  "corners", 
and  geometrical  optics  used  for  the  focussing  effects  of  curved  con¬ 
cavities. 

The  programs  described  in  this  paper  have  been  tested  and  have 
given  good  results  for  sectionally  continuous  targets  of  revolution 
having  straight  line  or  curved  smoothly-joined  profiles.  Extension 
of  these  programs  using  the  techniques  described  above  is  feasible. 
Thus  with  more  effort  a  computer  program  which  is  sufficiently 
general  to  treat  a  great  number  of  practical  targets  would  result. 


APPENDIX  I 

THE  GEOMETRIC  PROPERTIES  OF  A  GENERAL 
SECOND  ORDER  SURFACE  OF  REVOLUTION 

A  general  second  order  surface  may  be  described  in  spherical 
coordinates  by  the  following  equation 

(  103)  Ax  r2  sin20  cos2<j>  +  A2r2  sin20  sin2cj)  +  A3  r2  cos20 

+  A4  r2  sin20  cos  $  sin  c|>  +  A5  r2  sin  0  cos  0  sin  <j> 

+  A*  r2  sin  0  cos  0  cos  4>  +  Ayr  sin  0  cos  <)> 

+  A8r  sin  0  sin  <j>  +  A9  r  cos  0  +  Axo  r  sin  0  +  Axx  =  0 

For  a  surface  of  revolution  (he.,  a  surface  having  no  <j>  variation) 
we  can  reduce  Eq.  (  103)  to 

(  104)  Arl  r2sin20  +  Ar3  r2cos20  +  Bri  r2  sin  0cos  0 

+  Ar9  r  cos  0+  Arl0  r  sin  6  +  Aru  =  0 

by  examining  the  profile  of  the  surface  in  the  plane  4>  =  0.  For  the 
sake  of  convenience  Eq.(  104)  is  written  as 
(  105)  r2 U (  9)  +  r  V(  0)  +  ArU  =  0 


134 


135 


The  properties  of  this  surface  whicn  are  necessary  for  a  geo¬ 
metrical  theory  of  diffraction  analysis  are  the  unit  normal  vector 

A 

(  n)  ,  the  differential  arc  length  on  the  surface  (  ds)  ,  and  the  radii 
of  curvature  in  the  co-ordinate  directions  (  pj  in  the  0  direction,  p2 
in  the  (^-direction).  The  normal  vector  is  most  easily  obtained 
through  the  use  of  the  gradient.  Thus 


(  106) 

A  V [  r2U (  6 )  +  r  V(0)  +  Arll] 

n  =  - 

|v[  r2U (  0)  +  r  V(  0)  +  Arn  ]  | 

Taking  the  gradient  we  have 


(  107) 

r  (  ZrU  +  V)  +  0  (  rUe  +  V0) 

n  =  - T~ 

[  (  2r  U  +  V)  2  +  (  rU0+  Vq)2]2 

whe  re 

0U(  0) 

U0  =  00 

0  V(  0) 

0  00 

for  convenience  we  set 


(  108) 

Si  =  (  2r  U  +  V) 

S2  =  (  rU0+  V0) 

Now,  using  the  relations  between  unit  vectors  in  the  (  r,  8,40  and 


(  x,  y,  z)  systems. 


136 


A  /\  A  .  .  A 
(  109)  r  =  x  sin  0  cos  cj)  +  y  sin  0  sin  <\>  +  z  cos  0 

A  A  A  A 
0  =  x  cos  0  cos  <(>  +  y  cos  0  sin  §  -  z  sin  0 


we  may  write 


A  1  f  A  r  (  , 

(110)  n  =  t—  {x[Si  sin0cos<t>  +  S2cos0cos<t>J 


[sj  +  S*  ] 


+  y[Si  sin  0  sin  c|>  +  S2  cos  0  sin  <(>] 


+  z  [  Sj  cos  0  -  S2  sin  0]  } 


The  differential  arc  length  (  ds)  may  be  expressed  as  the 
scalar  product  of  the  derivative  of  the  position  vector  with  itself. 

For  a  curve  given  by  the  relation  f(  u,  v)  =0  we  have[  48] 

(111)  d F  =  Fu  du  +  Fv  dv 
and  the  differential  arc  length  is 

(112)  ds2  =  dF  •  d r 

This  relation  may  be  written  in  what  is  called  the  First  Fundamental 
Form  of  differential  geometry[  48]  viz. 

(113)  ds2  =  E  du2  +  2F  du  dv  +  G  dv2 


where 
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We  may  write  the  position  vector  in  rectangular  coordinates  as 


(114)  r  =  x  r  sin  0  cos  (|)  +  y  r  sin  0  sin  (()  +  z  r  cos  G 


The  0  and  <j>  derivatives  are  then 

(115)  rg  =  x[  rg  sin  0  cos  (j)  +  r  cos  0  cos  (()] 

+  y[  Tq  sin  0  sin  <j)  +  r  cos  0  sin  <\>] 

+  z[  rg  cos  0  -  r  sin  0] 

—  A  .  .  A 

(116)  r^  =  -x  r  sin  0  sin  4>  +  y  r  sin  0  cos  <j> 

The  derivative  of  r  with  respect  to  0  is  obtained  using  Eq.  (  105) 

dr(  6)  rZUe  +  rVe  rS2 

(117)  rfl  = -  =  -  -  -  -  - 

0  d0  2rU  +  V  Sj 

Next,  using  Eq.  (113)  and  algebraic  manipulation  one  obtains  in 
(r,0,<t>)  co-ordinates 

(118)  E  =  r2  +  rZ 
F  =  0 

G  =  r2  sin20 


as 


which  results  yield  the  arc  length  on  the  surface. 
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In  order  to  obtain  the  radii  of  curvature  on  the  surface  it  is 
necessary  to  evaluate  the  Second  Fundamental  Form  of  Differential 
Geometry[  48]  .  This  form  is  as  follows, 

(  119)  -dr  *  dn  =  e  du2  +  2f  du  dv  +  g  dv2 


in  this  equation  e,  f ,  g  can  be  computed  as 


e  = 

-  r 

u 

A 

"u 

f  = 

-7  • 
xu 

A 

"v 

g  = 

-  r 

A  V 

A 

nv 

Equation  (  119)  is  not  convenient  for  computation  in  the  case  being 
examined  as  it  generates  a  great  many  terms  which  must  be  reduced. 
An  equivalent  form  which  in  this  case  is  easy  to  evaluate  is  given  by 
Struik[  48]  as 


(  120) 


e  =  r 


uu 


f  =  7 


uv 


A 

n 


A 

n 


g  =  rvv  *  n 


This  form  reduces  the  number  of  terms  to  be  reduced  by  about  1/2. 


Using  Eq.  (115)  we  have 
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(121)  rgQ  =  x[r00sin0cos4>+2r0Cos0cos<i>-r  sin  0  cos  <j>] 

+  y[  x* g g  sin  0  sin  cj>  +  2rg  cos  0  sin  <j>  -  r  sin  0  sin  <j>] 
+  z[  r g g  cos  0  -  2rg  sin  0  -r  cos  0] 

_  A  A 


(  122) 

r04>  = 

-x  r  cos  0  sin  (j) 

+  y 

r  cos  0  cos  cj> 

(  123) 

A 

A 

r<M>  = 

-xrsinGcoscj) 

-  y 

r  sin  0  sin  (\> 

where 

►1 

CD 

CD 

II 

T0[  (  s  1  )  0  +  S2] 

+  r 

(Si)0 

Si 

(  s2  )  0  =  rgU  0  +  ruee+  vee 

Performing  the  operations  indicated  in  Eq.  (  120)  we  have,  after 
reduction 

Si  r00  +  2S2  r0  -  Si  r 

(  124)  e  = - 1 - 

[Si  +  s2]2 

f  =  0 

Si  r  sin20+  S2r  sin  0  cos  0 

g  = - t - 

[s{  +  s22]  2 

Now  the  normal  curvature  may  be  written  as[  41] 
e  du2  +  2f  du  dv  +  g  dv2 
E  du^  +  2F  du  dv  +  G  dv2 


(  125) 
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which  in  our  case  reduces  to 

e  d02  +  g  dcf>2 

'  # 

E  d02  +  G  d<t>2 

Now  the  curvatures  in  the  directions  of  the  coordinates  may  be  found 
by  taking  d0  and  d<J>  equal  to  zero  respectively. 

(  127)  Kj  =  e/E 

Kz  -  g  /G 

The  curvatures  Kj  and  K2  are  called  the  principal  curvatures.  The 
total  curvature  along  a  tangent  to  a  coordinate  direction  can  be 
expressed  by  introducing  the  angle  a  between  the  direction  d4>/d0  and 
the  curvature  direction  dcf>  =  0,  resulting  in 

(  128)  K  =  Ki  cos2a  +  K2  sin2a 

This  relation  which  expresses  the  normal  curvature  in  an  arbitrary 
direction  in  terms  of  the  principal  curvatures  is  known  as  Euler’s 


(  126)  Kn 


Theorem. 
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Two  further  expressions  are  useful.  These  are  expressions 
for  the  Mean  and  Gaussian  curvatures  and  are 
Eg  -2fF  +  eG 

(  129)  M  = -  ,  the  mean  curvature 

2(  EG  -  F2) 

and 

eg  -  f2 

(  130)  K  =  - 5-  ,  the  Gaussian  curvature. 

EG  -  F2 

As  the  radius  of  curvature  in  a  given  direction  is  the  reciprocal  of 
the  curvature  the  above  results  yield  the  radii  of  curvature  required 
in  the  evaluation  of  the  diffraction  and  attenuation  coefficients. 

The  above  results  allow  the  determination  of  the  radii  of 
curvature,  differential  arc  length,  and  Gaussian  curvature  at  any 
point  on  a  surface  of  revolution  defined  by  Eq.  (104).  These 
surfaces  include  the  figures  of  revolution  derived  from  the  conic 
sections  and  the  ogive. 
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APPENDIX  H 

DIFFRACTION  BY  A  PERFECTLY 
CONDUCTING  WEDGE 


In  order  to  satisfactorily  calculate  diffracted  rays  it  is  neces¬ 
sary  to  determine  the  diffraction  of  plane  and  cylindrical  waves  by  a 
wedge.  This  has  previously  been  done  (  see  Reference  26)  and  an 
expanded  treatment  is  presented  here.  Sommerfeld[  49]  obtained 
the  solution  for  a  perfectly  conducting  wedge  composed  of  two  half¬ 
planes  with  a  plane  wave  incident  on  the  wedge.  Sommerfeld  also 
obtained  an  explicit  form  of  the  solution  for  zero  wedge  angle;  i.  e.  , 
a  half-plane.  Pauli[  5]  determined  an  explicit  form  for  the  general 
wedge.  Oberhettinger[  13]  presents  a  different  form  for  the  general 
wedge  which  he  obtained  using  Green’s  function  techniques. 

The  diffraction  of  a  cylindrical  wave  by  a  wedge  is  obtained  by 
the  use  of  reciprocity,  together  with  Pauli’s  expressions  for  wedge 
diffraction.  Pauli's  expressions  give  the  diffraction  of  a  plane  wave 
by  a  wedge  for  a  general  angle  of  incidence  and  for  polarization 
either  perpendicular  or  parallel  to  the  edge  of  the  wedge.  The  total 
field  at  observation  point  P  of  cylindrical  co-ordinates  (  r,ijj)  ,  as 
shown  in  Fig.  50a,  is  given  by 
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INCIDENT 
CYLINDRICAL 
WAVE 

Xo-a> 


Fig.  50- -  Diffraction  by  a  wedge  of  included  angle  (  2-n)  tt. 

(131)  u(r,n,i(j)  =  v(  r ,  n,  Q)  -  v(  r ,  n,  iJj  +  \\>  Q) 

for  a  plane  wave  incident  from  direction  lJj  ,  The  plus  sign  applies 
for  the  polarization  of  the  electric  field  perpendicular  to  the  edge 


9u 

3n 


=  0 


wedge  f 

and  the  minus  sign  applies  for  polarization  parallel  to  the  edge 


(u 


wedg  e 


^  =0.  The  quantities  v(  r,  cj))  are  given  by 
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(  132) 


v.(  r ,  n;<j>)  =  v  (  r,  n,  <}>)  +  vg(  r ,  n,  4>) 


where  v‘  (  r,n,  4>)  is  the  geometrical  optics  field;  and  is  given  by 

f 


(13  3)  v'1'  (  r,  n,  <J>)  = 


where 


exp[  jkr  cos  ( <J>+  2irnN)  ]  , 

-TT+  2TrnN  <  4>  <  tt  +  2  it  nN 
0  otherwise 


N  =  0,  -  1  ,  -  2,  . 


vg(  r,n,4>)  is  the  diffracted  field  for  a  wedge  of  angle  (2-n)  it 
and  is  given  by 


V  akr 


where  a  =  1  +  cos  <j>. 


Equation  (  134)  is  composed  of  a  leading  term  plus  higher- 


order  terms  which  are  negligible  for  large  values  of  kr.  For  large 


values  of  (akr)  Eq.  (  134)  becomes, 


e-j(kr  +  W4)  I  sinZ 

(135)  VB(  r ,  n,  4>)  = - - - - ^ - - 

\l  2Trkr  cos  —  -  cos 

n  n 
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The  diffracted  field,  as  expressed  by  Eq.  (  135)  ,  is  that  from  which 
the  asymptotic  diffraction  coefficients  of  the  geometrical  theory  of 
diffraction  are  obtained[  l]  .  Since  this  expression  is  valid  only  for 
large  values  of  (  akr)  ,  it  is  not  valid  on  the  shadow  boundary  because 
a  =  0  there.  Then  Eq.  (  134)  must  be  used,  which  gives  the  value  of 
the  diffracted  field  on  the  shadow  boundary  as 

!  upper  sign  for  <j>  =  tt-* 

(  136)  vg(  r,  n,  tt)  =  +  y  e“J  r  +  •  •  •  {  #  + 

^  lower  sign  for  4>  =  tt 

The  value  of  v(  r,<j))  on  the  shadow  boundary  can  then  be  obtained 
from  Eqs.  (  132)  ,  (  133)  ,  and  (  136)  as 

(  137)  v(  r,n,  tt)  =  +  •  •  •  , 

which  is  one-half  of  the  incident  field  on  the  illuminated  side  of  the 
shadow  boundary. 

The  series  representation  of  vg  given  in  Eq.  (  134)  is  valid 
everywhere  except  for  the  values 

(138)  <t»  =  *  TT  +  2trnN,  n  4  2,  N  =  -  1,  ±  2,  ... 

Near  these  values  the  series  representation  converges  slowly,  and 
the  periodicity  property  of  the  exact  function  v(  r,n,cj>)  can  be  used 
to  overcome  convergence  difficulties  near  these  values.  The  exact 
function  v(  r,n,40  is  periodic  in  2irn  so  that 
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(139)  v(  r,  n,  4>  +  2imN)  =  v(r,n,4>),  N  =  0,  -  1,  -  2,  .  .  . 

Therefore,  if  the  series  representation  of  vg  converges  slowly 
because  cj>  is  near  one  of  the  values  expressed  by  N  4  0  in  Eq.  (  138)  , 
the  periodicity  property  of  Eq.  (  139)  can  be  used  to  represent  v(  r,<j>) 
by  employing  the  series  represenation  near  <j>  =  -  tt(  N  =  0).  The  only 
case  for  which  all  boundaries  are  regular  and  the  substitution  of 
Eq.  (  139)  is  not  necessary  is  for  the  thin  half-plane,  in  which  case 
n  =  2. 

We  now  examine  the  behavior  of  the  diffraction  coefficient 
vg(  r,n,4>)  for  the  case  illustrated  in  Fig.  51  in  which  two  reflected 
rays  exist  at  angles  ip  1  and  4^2*  We  first  examine  the  values  of 
4>  =  \\i  +  \\)  o  for  the  diffracted  rays  corresponding  to  the  directions  of 
these  reflected  rays.  The  values  of  the  pertinent  angles  are 

(  140)  4i  1  =  TT  -  4. o  ,  4*2  =  (  2nir  -  ir)  -  4»Q 

Thus  the  values  of  <t>  =  4J  +  40  are 

(  141)  4*  1  =  it  >  $2  =  (  ^niT  ”  w) 

For  ((>!  =  it  Eq.  (  136)  expresses  the  value  of  the  diffraction  coefficient 
and  thus  the  total  field  at  4  =  Is  one-half  the  geometrical  optics 
reflected  field  on  the  illuminated  side  of  the  shadow  boundary,  as 


given  by  Eq.  (  137). 
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Fig.  51--Case  of  diffraction  by  a  wedge  where 
two  reflected  rays  exist. 


However,  for  4>2  -  Znir  -  tt  and  n  4  2, 
and  consequently  the  following  substitution 


Eq.  (134)  is  not  analytic 
is  made: 


(142)  v(  r ,  n,  <)>)  “•*  v(  r,  n, 

for  values  of 


(  143) 


>  Trn 
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The  region  over  which  the  substitution  of  Eq.  (  142)  is  made,  as 
expressed  by  Eq.  (  143)  is  determined  from  the  symmetry  property 
of  v(r,n,  ()));  i  e.  , 

(144)  v(r,n,  -<j>)  =v(r,n, +<|>) 

That  is,  if  the  substitution  of  Eq.  (  142)  is  made  for  values  of  cj)  <  im, 
then  the  substituted  values  of  -  2im  are  closer  to  the  value  of  non- 
analyticity  Tr-2mr  than  the  original  values  of  were  to  the  value  of 
nonanalyticity  2mr-Tr.  It  is  also  noted  that  the  value  substituted  for 
v(  r ,  n,  <j>2)  is  v(  r,  n,  -t r)  ,  which  corresponds  to  the  correct  value  for 
a  shadow  boundary. 

The  equations  given  in  this  Appendix  apply  for  plane-wave 
incidence,  but  they  can  be  used  for  cylindrical-wave  incidence  as 
shown  below.  Consider  the  two  situations  shown  in  Fig.  50.  It  is 
desired  to  find  the  field  ua  in  some  direction  £  for  the  wedge  illumi¬ 
nated  by  a  cylindrical  wave  with  its  source  located  at  (x0,  or).  By 
reciprocity,  the  field  ua  is  equal  to  the  field  u^  which  is  located  at 
the  point  (  r  =  xQ,  i[>  =  a  )  and  with  a  plane  wave  incident  from  the 
direction  =  The  value  of  u^  is  given  by  Eq*  (  134).  Thus 
using  the  property  expressed  in  Eq.  (  131)  the  solution  for  dif¬ 
fraction  of  a  cylindrical  wave  by  a  wedge  becomes 


(  145) 


ua  =  v(  xo>  n>  £+  )  -  v(  x0,  n  |  -  a  ) 
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The  field  at  infinity  is  given  by  Eq.  (145)  for  a  perfectly  conducting 
wedge  illuminated  by  a  line  source  at  (xQ,a). 

If  the  quantity  (akxQ)  is  sufficiently  large  for  Eq.  (135)  to  be 
valid,  then  the  diffraction  pattern  for  cylindrical- wave  incidence  has 
the  same  form  as  that  for  plane-wave  incidence.  In  other  words, 
the  diffraction  for  cylindrical- wave  incidence  is  the  same  as  that 
for  plane-wave  incidence  in  regions  sufficiently  removed  from  the 
shadow  boundaries.  The  region  near  the  shadow  boundary  in  which 
Eq.  (134)  must  be  used  increases  for  decreasing  values  of  xQ  and 
may  encompass  all  360  degrees. 

Hutchins  [  5 0  ]  has  generalized  the  Pauli  solution  presented 
above  and  has  obtained  a  series  solution  which  is  valid  for  the 
exterior  and  interior  regions  of  a  wedge.  In  addition  his  solution 
is  accurate  for  wedges  of  large  included  angles  where  the  Pauli 


solution  is  not  accurate. 
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APPENDIX  HI 

APPROXIMATION  FOR  THE  SPECULAR  SCATTERING 
BY  A  GENERATOR  OF  A  CONE 


Specular  scattering  by  a  generator  or  frustrum  of  a  cone  is 
shown  in  Fig.  52.  An  approximate  solution  for  the  backscattered  field 


Fig.  52--Scattering  by  a  conical  frustrum. 

in  and  near  the  specular  direction  may  be  constructed  as  follows. 

It  is  assumed  that  thebackscattered  field  close  to  the  normal  to  the 


cone  can  be  written  as 
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(146)  US(U  =  +±i_  F(r)  e-Jkr  .  ^ 

\J2u  k  kd  sin  (  £  )  r 

Equation  (146)  combines  the  sin(x)/x  pattern  function  behavior  previ¬ 
ously  discussed  and  the  spatial  attenuation  factor  F{r).  Proceeding 
from  Fig.  5Z  we  have  that 

( 147)  px  sc  ai  / cos  rj 

(148)  =  a2/ros  r) 

where  pl  and  are  the  finite  radii  of  curvature  at  the  ends  of  the 
frustrum. 

11491  F(r)  U=o'Ji 


where 

(150) 


P  =  Pi  +  T~  (p2  "  pi>  • 


The  average  spatial  attenuation  for  £  small  is  approximated  as 

(151)  F(r)  ave  =  4-  \  \l - 4 -  dz 


-  i-  D 


with  the  result 


(162)  F(r ) 


ave 


n/2  1 

3  P2  -  Pi 


(p23/2-  Pl3/i]x  . 


The  approximate  scattered  field  from  the  generator  of  the  cone  is  then 


(153) 


USapp  (?)  =  +  -“T - CMase  3i.1kd  T  0  ^  •  U1 


\[2u  k 


kd  sin  £ 


APPENDIX  IV 

THE  DETERMINATION  OF  THE  GEODESIC  PATHS 
ON  A  GENERAL  QUADRIC  OF  REVOLUTION 


The  determination  of  the  geodesic  paths  on  a  general  quadric 
surface  of  revolution  is  not  an  easy  task  in  general.  It  is  known 
that  the  geodesic  paths  on  a  body  of  revolution  can  be  found  by 
quadratures.  However,  in  the  case  studied  here  it  is  convenient 
to  utilize  the  differential  equation  for  the  geodesic  given  by 
Struik[  48]  as 

(154)  v"=A(vi)3+B(v,)2  +  cvi  +D 

G0 

where  A  =  -  - 

2E 

Ee  Ge 

B  “  E  "  2G 

Ee  Ge 

2E  G 


D  = 

2G 


v  =  <|>(e) 

and  where  E  and  G  are  given  by  Eq.  (118)  and  Eq  and  Gq  are  the 
partial  derivatives  with  respect  to  0. 
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We  now  assume  a  Taylor  series  solution  for  v  =  <j)  in  the  form 


The  differential  equation  is  then 


And  the  following  procedure  will  yield  a  numerical  solution  by 
interation: 

1.  Determine  d<^»  /dG | 

ei 

2.  Compute  A,  B,  C,  D,  at  Gp  <\>i 

3.  Compute  d2< 4>  /d02  I  .  n  , 

( «t»i) 

4.  Compute  <\>(  )  using  the  series 

5.,  Compute  d<j>/dG  =  (  -  <{>•)  /  A0 

6„  i  =  i+1  and  proceed. 
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APPENDIX  V 

THE  COMPUTER  SOLUTION  FOR  BACKSCATTER 
BY  A  SEC TIONALLY  CONTINUOUS  BODY  OF 
REVOLUTION  DESCRIBED  BY  A  SECOND 
DEGREE  EQUATION 

In  the  computer  solution  Eq.  (  104)  is  used  to  describe  the 
profile  of  the  target  within  a  section  defined  by  the  angles  THTI(  I) 
and  THTF(I).  The  constants  must  be  arranged  so  that  Eq,  (  110) 
generates  an  outward  normal  to  the  surface.  The  angle  of  the 
incident  wave  is  THT  and  the  wavelength  is  WAVE.  The  following 
flow  diagrams  describe  the  computer  programs  which  were  coded 
in  Fortran  IV  for  the  solution  of  the  backscattered  fields  of  the 
target.  A  list  of  the  Subroutines  and  Functions,  other  than  standard 
library  routines  which  were  used  is  also  presented. 

The  programs,  titled  "Wedge  Diffraction  Computer  Program" 
and  "Creeping  Wave  Computer  Program"  utilize  the  analytic  tech¬ 
niques  presented  previously.  The  Wedge  Diffraction  Computer 
Program  uses  the  wedge  diffraction  coefficient,  the  diffraction  by  a 
pair  of  wedges,  and  the  solution  for  the  diffracted  fields  near  the 
axis  to  construct  the  solution  for  the  scattered  fields  due  to  single 
diffraction  from  locations  on  the  target  where  the  slope  is  discontinuous. 


The  creeping  wave  computer  program  uses  the  geometrical  optics 
field  from  the  specular  point  and  the  creeping  wave  solution  for 
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waves  which  propagate  around  the  target  to  construct  the  solution 
for  a  target  having  no  discontinuities  in  slope  on  its  surface.  The 
scattered  fields  for  a  target  where  both  wedge  diffraction  and 
creeping  wave  contributions  exist  can  be  constructed  using  a 
combination  of  the  results  of  both  programs. 
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Wedge  Djffract/o n  Computer  Program 
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SUBROUTINES  AND  FUNCTIONS 

FCOMM(  I)  Shifts  constants  describing  each  section  from  /DATA/ 
COMMON  to  unlabeled  COMMON,  thus  adjusting  the  constants 
to  the  current  section  (  I). 

RAD(  THT)  Computes  the  radial  distance  to  the  surface  at  the  angle 
THT. 

FNORM(  FNVX ,  FNVY,  FNVZ,  R,  THT,  PHI)  Computes  the  surface 
normal  N  =  xFNVX  +  yFNVY  +  zFNVZ  at  (  R,  THT,  PHI) 

WILLY  (  FWI,  THTI,  PHI,  THIN,  PH  IN)  Indicates  if  the  location 

(  RAD(  THTI)  ,  THTI,  PHI)  is  illuminated  (  FWI  =  1.  )  or 
shadowed  (FWI-0.  )  with  respect  to  the  incident  wave 
(  THIN,  PHIN). 

CROSS(  X,  Y,  Z,  Al,  A2,  A3,  Bl,  B2,  B3)  Computes  X  -  AXB. 

PHASE(  THTI,  PHII,  THTB,  RB,  FK)  Computes  the  backscattered 
phase  of  the  location  (  RB,  THTB)  for  a  wavenumber  (  FK) 
and  incident  wave  (  THTI,  PHII). 

SINXX(  Y)  Computes  sin  (  Y)  /Y  . 

sin  (  tt/fN) 

DIFF1(  FN,  PHI,  BETA)  Computes  DIFF1  = - 

FN*  sin(  BETA) 


c  os(  tt/FN)  -  COS  (  phi/fn) 
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BESL1(  X)  Computes  Jj  (  X) . 

FSPDT(  RSP,  THSP,  FNT,  THTS,  THTF,  DTHT,  PHI,  VX,  VY,  VZ) 
Increments  by  steps  of  DTHT  within  the  limits  THTS  to 
THTF  and  finds  the  location  (  RSP,  THSP)  at  which  the  scalar 
product  FNT  =  N  •  V  is  a  maximum. 

FINT(  SSS,  FCTI  ,FLL,  FUL,  ERRR,  NX)  Performs  a  numerical 
integration  of  the  complex  function  FCTI  between  the  limits 
FLL  to  FUL,  within  a  percent  error  ERRR,  returning  the 
answer  in  SSS.  If  NX  is  1  equal  increments  are  used,  if  2 
adjusted  increments  are  used. 

DELSP(  THT)  Computes  the  incremental  arc  length  at  the  angle 
THT. 

SPHASE(  THTI,  PHII,  THTB,  RB,  FK)  Computes  the  incident  phase 
of  the  location  (  RB,  THTB)  for  a  wavenumber  (  FK)  and 
incident  wave  (  THTI,  PHII) . 

ALPH(  RHOl,  RH02,  WAVE)  Computes  the  attenuation  function 
(ALPH)  for  radii  of  curvature  (  RHOl,  RH02)  and  a  wave¬ 
length  (  WAVE)  . 

DSQ(  RI1,  RI1,  WAVE)  Computes  the  square  of  the  diffraction 
coefficient  (  DSQ)  for  radii  of  curvature)  RI1,  RI2)  and  a 


wavelength  (WAVE). 
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ALPHDS(  THT)  Computes  the  product  of  the  attenuation  coefficient 
(  ALPH)  and  the  incremental  arc  length  (  DELSP)  at  the  angle 
(  THT). 

DIFFGO(  R,  THT,  RTHT,  RTTH,  ECAP,  FCAP,  GCAP,  ELC,  FLC,  GLC) 
Computes  the  radial  first  (  RTHT)  and  second  (  RTTH) 
derivatives  of  the  surface  and  the  differential  geometry 
parameters  (  ECAP,  FCAP,  GCAP,  ELC  ,  FLC  ,  GLC)  at  a  point 
on  the  surface  (  R,  THT). 
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